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ABSTRACT 

In this third paper in a series of three, we present a detailed study of the broadband 
spectral energy distribution (SED) of active galactic nuclei (AGN) based on a nearby 
unobscured Type 1 AGN sample. We perform a systematic cross-correlation study of 
several key parameters, i.e. T2-iokeV, i2-iofceV, Ltoi/LEdd = ^Edd, Lboi/L2-wkeV = 
K2^iokeV, Lboi/LbwoA = i^zwoA-, FWHM//^ , Mbh, ctox-, ax and auv ■ The weh de- 
fined spectral properties of the sample enable us to improve existing relations and to 
identify new correlations among these parameters. We confirm a break region around 
FWHM/f/3 ~ 4000 km s'^ in the V2-wkeV vs. FWHMff/g correlation and hog(MBH) 
~ 8.0 in the T2-wkeV vs. y^BH correlation, where these correlations appear to change 
form. Beyond the break point the intrinsic T2~wkeV index is dispersed around 1.8. 
Several new correlations are also reported in this paper e.g. strong correlations in 
K5100 vs. \Edd, K5100 vs. K2^wkeV and K2-iofeey vs. Mbh- The principal component 
analysis (PCA) is performed on the correlation matrix of the above parameters. This 
shows that the three physical parameters, i.e. black hole mass, mass accretion rate 
and Eddington ratio, drive the majority of the correlations. This is consistent with 
PCA results found from previous optical spectral studies. 

For each key parameter, we split the AGN into three sub-samples, binned based 
on increasing value of that parameter. We co-add the model SEDs for each object in 
the sub-sample to see how the SED changes with that parameter. Most parameters, 
except Ijboh show similar systematic changes in the SED such that the temperature 
at which the disc peaks is correlated with the ratio of power in the disc versus the 
Comptonised components and the hard X-ray spectral index. This underlying change 
in SED shape shows that AGN do exhibit intrinsically different spectral states. This 
is superficially similar to the SED differences in BHB (black hole binary) seen as A Edd 
increases, but the analogy does not hold in detail. Only objects with the highest \Edd 
appear to correspond to a BHB spectral state (the disc dominated high/soft state). 
The AGN with typical mass accretion rates have spectra which do not match well with 
any state observed in BHB. We speculate that this could be due to the presence of 
a powerful UV line driven disc wind, which complicates simple mass scaling between 
stellar and supermassive black holes. 
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1 INTRODUCTION 

The broadband spectral energy distribution (SED) of ac- 
tive galactic nuclei (AGN) has been studied for many years. 
The unobservable gap between the extreme UV and soft 
X-rays, imposed by Galactic photoelectric absorption sets 
a major barrier to reconstructing the overall SED. This is 
particularly problematic because it is in this energy range 
that the peak of the SED often occurs. Due to this diffi- 
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culty, most previous work focused on spectral properties of 
either the optical/UV side (e.g. Sanders et al. 1989; Francis 
et al. 1991; Zheng et al. 1997; Vanden Berk et al. 2001), 
or the X-ray side (e.g. Wilkes & Elvis 1987; Green et al. 
1995; George et al. 2000). However, the multi- wavelength 
combination is required to see the overall SED behaviour 
(e.g. Elvis et al. 1994; Laor et al. 1997; Richards et al. 2006; 
Shang et al. 2005; Lusso et al. 2010, hereafter: LussolO). 
In order to recover the spectral shape in the UV/soft X- 
ray gap, some studies simply connected the optical/UV and 
X-ray spectra to estimate the SED (e.g. Elvis et al. 1994; 
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Richards et al. 2006; Shang et al. 2005; LussolO). As the ori- 
gin of the optical and hard X-ray emission became clearer, 
more broadband SED models with multiple physical com- 
ponents were used to fit the multi-waveband spectra. These 
were used to extrapolate across the UV gap based on the 
model assumptions, but with better justification than be- 
fore (e.g. Middleton, Done & Gierlihski 2007; Vasudevan & 
Fabian 2007, hereafter: VF07; Vasudevan & Fabian 2009, 
hereafter: VF09; Jin et al. 2009; Grupe et al. 2010, here- 
after: GrupelO). For example, VF07 used a power law plus 
accretion disc model to construct the broadband SED for 
their sample. This technique was repeated in VF09 based 
on 29 AGNs from the reverberation mapping sample (Pe- 
terson et al. 2004), and using simultaneous optical/UV and 
X-ray data from XMM-Newton. GrupelO used an absorbed 
power law for the X-ray spectrum, and an exponentially cut- 
off power law (or a broken power law) for the optical/UV 
spectrum, to construct the broadband SED for their sample, 
for which simultaneous optical, UV, X-ray data are avail- 
able from SWIFT. Further discussion of some previous SED 
studies is given in Jin et al. (2011) (hereafter: Paper-I). 

Within these works, some specific AGN parameters 
were studied, including black hole mass (Ms//), Edding- 
ton ratio {L/LEdd = ^Edd), bolometric luminosity (Ltoi), 
2-10 keV luminosity (L2-iofcev), 2-10 keV photon index 
(r2-iofcev), H/3 FWHM. Also, some parameters directly re- 
lated to the broadband SED shape were proposed, such as 
the optical to X-ray spectral index (qoJ3, Tananbaum et al. 
1979), 2-10 keV bolometric correction (K2-iofceV, defined as 
L6oi/L2-iofceV, e.g. VF07), 5100A luminosity scaling factor 
(ksioo, defined as Lboi/Lsioo, where L5100 is the monochro- 
matic luminosity at 5100A, Kaspi et al. 2000). 

Since the number of AGN with both high quality op- 
tical/UV and X-ray spectra is relatively small, much effort 
is devoted to search for correlations among the key SED 
parameters, especially those parameters capable of repre- 
senting the broadband SED shape. Then, for those many 
AGNs lacking sufficient spectral information, these correla- 
tions can be used to predict the SEDs that cannot be defined 
from direct observation. Indeed, many such parameter cor- 
relations have been proposed. For example, VF07 reported a 
strong correlation between «:2_iofceV and Xsdd (VF07; VF09; 
LussolO). Correlations were also found between H/3 FWHM 
and r2-iofceV (e.g. Leighly 1999; Reeves & Turner 2000; 
Shemmer et al. 2006, hereafter: S06; Shemmer et al. 2008, 
hereafter: SOS; Zhou & Zhang 2010, hereafter: ZhoulOa), be- 
tween Xsdd and r2-iofceV (e.g. Lu & Yu 1999; Porquet et al. 
2004; Wang, Watarai & Mineshige 2004; Bian 2005; S06,08), 
and between H/3 FWHM and Xsdd (e.g. GrupelO; Paper-I). 
The correlation between L2500 and L2fcev is another impor- 
tant result, which led to further correlation studies related 
to Oox, L2500, L2fceV and redshift (e.g. Green et al. 2009, see 
references given in Table [2^ . In table lAll we list some of the 
principal parameters, and list the relevant papers discussing 
the correlations. Later in the papers we will discuss some 
additional correlations e.g. Mbh vs. T2-iokeV and Mbh vs. 
K2-iofceV, in the context that changes in these parameters 
are caused by changes in the fundamental physical processes. 

^ ciox = — ^"^^^^'^^"^^^^^^^^^^^^ ^ assuming F(i/) oc to en- 
sure a non- negative aox- 



However, due to the difficulty in obtaining and ana- 
lyzing both optical/UV and X-ray spectra for a large sam- 
ple, these parameter correlations are reported separately and 
based on different samples, rather than being studied sys- 
tematically for a single well-defined sample. Furthermore, 
most of the previous samples were not selected based on 
their spectral properties, so effects such as reddening will 
introduce biases into the cross-correlations. The lack of a 
self-consistent physically motivated broadband model, has 
also been a constraint to performing a more detailed SED 
study. 

In Paper-I, we defined a sample of 51 Type 1 AGNs 
with both optical/UV and X-ray spectra which are of high 
quality, and without evidence of complex spectral absorption 
features e.g. a warm absorber. Based on this bright and un- 
obscured Type 1 AGN sample, we applied our latest optical 
and broadband SED model to perform the spectral fitting, 
and so matched the optical spectrum and produced a broad- 
band SED for each AGN in the sample. This is so far the 
most detailed spectral analysis for a medium sized sample 
of AGNs, with such well defined high quality spectra. In our 
preceding paper (Jin, Ward & Done 2012, hereafter: Paper- 
II), we studied the Balmer emission line properties and the 
Relation between the optical and X-ray emission, based on 
the same sample and the spectral fitting results given in 
Paper-I. This paper is based on the sample of Paper-I, but 
focuses on the shape of the broadband SED, and how this 
relates to a range of parameters from the model. We will 
approach this problem by first investigating the numerous 
correlations previously reported. Then a set of mean SEDs 
based on key parameters are constructed and studied in de- 
tail. 

These mean SEDs show a clear change in shape as a 
function of the AGN parameters, probably most fundamen- 
tally driven by a change in Xsdd- The SED becomes more 
disc dominated and the X-ray tail softens as Xsdd increases 
(e.g. Grupe et al. 2010). This is superficially similar to the 
changes seen in the Galactic black hole binary (BHB) sys- 
tems as the mass accretion rate increases. In BHB at low 
Eddington ratio (A^dd =^ 0.02) the X-ray spectra are hard, 
< 2, and the disc emission is weak (low/hard 
state). As the mass accretion rate increases the disc increases 
in importance relative to the tail, and the tail steepens until 
it reaches the disc dominated 'thermal state' (also known as 
the high /soft state) where the power law tail is somewhat 
steeper r2_iofceV ~ 2 — 2.2. At high Eddington fractions the 
source can also show a 'very high' or steep power law state, 
where both disc and tail contribute to the spectra, but the 
spectra are typically steep (r2-iofceV 2.2) (Remillard & 
McClintock 2006; Meyer-Hofmeister, Liu & Meyer 2009). It 
has been proposed that AGNs are simply scaled up counter- 
parts to the BHBs, so they should also show similar spec- 
tral states (Done & Gierliiiski 2005; McHardy et al. 2006; 
Gierlihski et al. 2008; Middleton et al. 2009; Jin et al. 2009). 
Therefore the Eddington ratio may also be a good indica- 
tor an AGN's accretion state, and so determine the shape of 
its SED, explaining the multiple strong correlations between 
various SED parameters. 

This paper is organized as follows. Section 2 gives a 
brief description of the main characteristics of the sample, 
and the method used for fitting the broadband SED. Fur- 
ther details of this are presented in Paper-I. In Section 3 
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Table 1. The mean parameter values with one standard deviation for our sample, together with some recently published samples. N: 
sample size. K2-10: the 2-10 keV bolometric correction. LussolO: Lusso et al. (2010); GrupelO; Grupe et al. (2010); ZhoulOa; Zhou & 
Zhang (2010); VF07,09; Vasudcvan & Fabian (2007, 2009). m: the Eddington ratios were calculated using the luminosity dependent 
2-10 keV bolometric correction in Marconi et al. (2004); r: the black hole masses are from the reverberation mapping study in Peterson 
et al. (2004). *: there are 114 sources in ZhoulOa, but we only consider the 69 NLSls and BLSls with both FWHM and \Edd 
measurements. 



N Redshift T2-iok&v K2-10 ksioo >^Edd FWHMif^j Mbh o-ox ^boi 

Sample km log(MQ) log{erg s~^) 



This Work 


51 


n 1 Q7-t-0. 158 


1.91±0.26 




I5tl* 


27+°-^l 


Qc;Kn+3880 
■JJ"U_lg60 


7.93±0.52 


1.35±0.14 


45.47±0.57 


LussolO 


545 






97-H28 
^'-14 










1.40±0.16 


45.54±0.57 


GrupelO 


92 


0.112±0.077 






21±3 


1.87±3.26 


2670±1670 


7.37±0.73 


1.42±0.17 


45.00±0.96 


ZhoulOa 


*69 


n 050+°-l°3 

U.UOU_g 


1.97±0.29 






mn 9/1 +0.76 
'^■^^-0.18 


2600+2500 
/:duu_^280 








VF09 


29 


„ „oo+0.074 
U.UOO_Q 


1.85±0.32 






0.18±0.16 




'•7.93±0.66 


1.39±0.24 


44.89±1.00 


VF07 


54 


064+° "'^ 




2g+39 
^"-16 




U-^0_0.13 




7.89±0.82 




45.20±1.01 



we examine the parameter correlations across five parame- 
ter groups. Some significant new correlations are proposed. 
In Section 4 we perform a systematic study of all correla- 
tions among several selected key parameters, and build a 
cross-correlation matrix. The principal component analysis 
technique is used on this correlation matrix in order to derive 
the eigenvectors. In Section 5 various mean SEDs are con- 
structed, based on the mean values of the key parameters. 
In Section 6 we investigate further the reliability of these 
correlations, including the effect of a correction to the black 
hole mass from radiation pressure. In Section 7 we summa- 
rize our results, and propose topics for further study. A flat 
universe model with a Hubble constant of Ho = 72 km s~^ 
Mpc~i, 0,M = 0.27 and CIa = 0.73, is adopted throughout 
the paper. 



2 PROPERTIES OF THE SAMPLE AND 
BROADBAND SED MODELING 

2.1 Sample Selection and Properties 

The AGN sample used in this paper is reported in Paper-I, 
based on the cross-correlation of SDSS DR7 and XMM- 
Newton catalogs. The complete source list, detailed sample 
selection criteria and sample properties can be found in 
Papcr-I. The main selection steps are summarized below: 

(1) We searched the 2XMMi and SDSS DR7 catalogs and 
identified 3342 extragalactic sources having both X-ray and 
optical spectra. 

(2) Within these sources, wc selected those with HjS in 
emission and redshift z < 0.4, so that both the Ha and 

emission lines are covered by the SDSS spectra. This 
assists with modelling of the Balmer lines (see Paper-I). 
This selection resulted in 802 unique X-ray sources. 

(3) Within this subsample, we identified 96 Type 1 AGNs 
all with a minimum of 2000 counts in at least one of the 
three XMM-Newton EPIC cameras, to ensure high X-ray 
spectral quality. 

(4) Wc then excluded 23 sources whose Hfi line was 
modified due to high reddening, low S/N or a data gap in 
the SDSS spectra. The resulting sample contains 73 AGNs. 

(5) For each of the 73 sources, a power law model was 
fitted to the rest-frame 2-10 keV X-ray spectra. The 16 
objects with photon index uncertainties greater than 0.5 



were excluded, leaving 57 Type 1 AGNs with relatively well 
constrained 2-10 keV spectra. 

(6) A further 6 objects were excluded because of an 
obvious signature of a warm absorber at ~0.7 keV. 
This criterion means that the observed spectra are very 
likely directly related to emission from the bare central core. 

The final sample contains 51 AGNs, with 12 AGNs clas- 
sified as NLSl using the conventional definition (Goodrich 
1989), while the others are all BLSls. The BAL quasar PG 
1004-1-130 is also included in our sample whose weak and 
featureless X-ray emission is still under debate (Miller et al. 
2006) The vast majority of the sample are radio-quiet with 
only three radio-loud sources i.e. PG 1004-1-130, RBS 0875 
and PG 1512-^370. 

The FWHM of the sample ranges from 600 km s"i to 
13000 km s-^ 



2.2 The Spectral Modelling 

Based on the high quality spectra, Paper-I conducted de- 
tailed spectral analysis for each source in the sample. In 
the optical spectral fitting, the Ha and H/3 lines were fit- 
ted using three components. The narrow component has 
the same jjrofile as the entire [OIII]A5007 line, i.e. includ- 
ing both the central and blue component in [OIIIJA5007. 
The intermediate and broad components are assumed have 
the shape of a Gaussian. All other strong nearby emission 
lines are included by adding more Gaussian profiles into the 
whole model. Then a complete model with multiple com- 
ponents were used to fit the whole SDSS spectra, including 
the underlying continuum approximated by a power law, the 
Balmer continuum, the Fell 'false' continuum and all strong 
emission lines. The broadband SED model used in Paper-I 
{optxagn, hereafter: Model- A) consists of the following three 
continuum components: 1. emission from a modified stan- 
dard accretion disc, whose energy within the corona radius 
is completely reprocessed into the other two high energy 
Comptonisation components. 2. emission from the low tem- 
perature, optically thick Comptonisation, which mainly ac- 
counts for the soft X-ray excess. 3. emission from the high 
temperature, optically thin Comptonisation which gives the 
power law shape of the hard X-ray spectrum above 2 keV. 
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SED Example 
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Figure 1. An example of the broadband SED fitting using op- 
txagnf model in Xspec vl2 which includes the effect of a colour 
correction. The data is taken from SDSS and XMM-Newton ob- 
servations of RBS 769. The solid red line shows the total model; 
the dashed green line shows the colour corrected and truncated 
accretion disc emission; the dotted orange line shows the low tem- 
perature optically thick Compotonisation; the dot-dash blue line 
shows the high temperature optically thin Compotonisation. The 
reduced is 1.16 for this spectral fitting. 



Both Galactic extinction and the small amount of intrinsic 
reddening/extinction are included in the model. 

Because the model was in the process of development, 
Paper-I did not consider the effect of a colour temperature 
correction in the accretion disc model (e.g. Ross, Fabian & 
Mineshige 1992; Davis & Hubeny 2006). This effect is due to 
the fact that the absorption opacity decreases significantly 

as the black hole IHSjSS illCr63,S6S (/^abs ^ 

M"^/*, Done et 

al. 2011), so for the same Eddington ratio, AGN's disc has 
a lower absorption opacity than BHB. Then for some AGNs 
the accretion disc may no longer be locally thermalised. The 
higher temperature photons can emerge from regions deeper 
in the disc, and so the disc spectrum extends to higher en- 
ergy than for standard accretion disc spectrum, producing 
the effect of a colour temperature correction. However, this 
can only happen for disc regions ofT > 3x10* K where suf- 
ficient hydrogen atoms are ionized. In addition, the electron 
scattering opacity also becomes important in these regions. 
The lower energy electrons in the surface regions Compton 
down-scatter the higher temperature photons from deeper 
disc regions, thus reducing the colour temperature correc- 
tion. These two competing effects lead to an effective colour 
temperature correction (Davis et al. 2006). 

The maximum effective temperature of the accretion 
disc is kT ~ 10(m/M8)^/'' eV (where m = Ltoi/LEdd, 
Ms = M/10* Mq), so only for AGNs with low mass black 
holes and high mass accretion rates such as the NLSls (e.g. 
Boiler, Brandt & Fink 1996; Paper-I), electron scattering 
opacity will dominate, and the effect of colour tempera- 
ture correction becomes important. A typical colour tem- 



perature correction of 2.6 is predicted for an AGN with 
Mbh^IQ'^Mq, Ai5dd=l-0 (Davis, Done & Blaes 2006; Done 
et al. 2011). This effect combines with an already hot disc 
due to the low black hole mass, resulting in a disc spectrum 
that extends significantly into the soft X-ray range. 

A more detailed explanation of this colour temperature 
correction in AGN can be found in Done et al. (2011), where 
a model {optxagnf in Xspec vl2) is proposed which applies 
a colour temperature correction for the Comptonised accre- 
tion disc model (hereafter: Model-B). In this paper we use 
this more advanced model to fit the optical, UV and X-ray 
data so as to construct the broadband SED for each source. 
Figure [1] shows an example of the broadband SED fitting 
to the multi-waveband spectra of RBS 769. A detailed de- 
scription of Model-B in comparison with the earlier model 
optxagn (Model- A), can be found in Done et al. (2011). The 
set of fitting parameters for each source, based on Model-B, 
is given in Table [Bl] Section [6.31 discusses the statistical dif- 
ferences resulting from using Model-B and Model-A for our 
sample. 

All the principal SED parameters such as XEdd, 
1^2-iokeV and Oox are calculated from the new model fit- 
ting (see Table [B2| . A cross-correlation study of the various 
parameters is then conducted. In Table [T] we compare the 
mean values of some SED parameters of our sample with 
those samples used in previous work. The result of correla- 
tions established in these previous studies will be compared 
with ours in the following sections. 

There are two sources that we treated as being anoma- 
lous in our study. The first is PG 1004-f 130, a broad absorp- 
tion line (BAL) quasar, whose X-ray was reported as being 
extraordinarily weak. Although its X-ray spectrum does not 
show clear absorption edges, it is nevertheless likely to be 
heavily absorbed, so that the remaining X-rays may have a 
different origin such as a sub-parsec-scale jet. It has been 
suggested that the X-ray emission from PG 1004+130, af- 
ter correcting for its intrinsic absorption, is 0.73 dex weaker 
than normal PG radio loud quasars (PG RLQs) normalized 
to similar optical/UV luminosities (Miller et al. 2006). Due 
to its distinct X-ray spectrum (and correspondingly differ- 
ent X-ray parameters), we did not include this source in our 
regression analysis. The other anomalous source is Mrk 110. 
This source shows strong optical variability (KoUatschny et 
al. 2001; KoUatschny 2003) and its SDSS spectrum has a 
different slope from the (non-simultaneous) XMM-Newton 
OM data (see Paper-I). However, the optical spectrum ob- 
tained using FAST shown in Landt et al. (2011), is consistent 
with the OM data, and is also an order of magnitude more 
brighter than the SDSS spectrum. Therefore the SDSS spec- 
trum of Mrk 110 is not consistent with its broadband SED 
parameters, and so we exclude Mrk 110 from any correlation 
which depends on optical parameters, but include it for any 
UV/X-ray correlations. 



3 INVESTIGATING THE CORRELATIONS 
FOR VARIOUS SED PARAMETER GROUPS 

In this section we divide the SED parameters into several 
sub-groups, based on correlations reported in the literature 
(see Section[T]), and then perform a cross-correlation analysis 
within each group. 
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Table 2. Comparison of regression line coefficients for L2fcey, L250O1 <^ox and redshift Correlations. LussolO: Lusso et al. (2010); 
GrupelO; Grupe et al. (2010); GreenOQ: Green et al. (2009); Just07: Just et al. (2007); SteffenOe: Steffen et al. (2006); Strateva05: 
Strateva et al. (2005); HasingerOS: Hasinger (2005). opt: optically selected sample, xray: X-ray selected sample, ^linear regression 
results using lj2ke.v L2500 from the reconstructed broadband SED corrected for both intrinsic and Galactic reddening/extinction, 
^linear regression when lj2keV ^'^d L2500 were not corrected for the best-fit intrinsic reddening/extinction, but corrected for the Galactic 
value, "for the SDSS main sample -|- high-z sample -|- Sy 1 Sample (see Strateva05). * measured directly from the regression line in 
Fig. 5(b) of Hasingcr05. 
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3.1 Group 1: L2500J L2fceV and Uox 

The aox index has been used as the indicator of the SED 
shape for more than thirty years. It is often used as an 
indication of the AGN's broadband SED, and to convert 
between the AGNs' optical luminosity function (OLF) and 
X-ray luminosity function (XLF) (e.g. Hopkins, Richards & 
Hernquist 2007, hereafter: Hopkins07 ; Tang, Zhang & Hop- 
kins 2007, hereafter: Tang07). 

Many studies have been carried out on the evolution of 
aox with luminosity (e.g. Avni & Tananbaum 1982; Wilkes 
et al. 1994; Strateva et al. 2005; Steffen et al. 2006; Just 
et al. 2007; Green et al. 2009, hereafter: Green09; LussolO; 
GrupelO), which may provide clues on the emission mecha- 
nism. The value of aox has been found in the range 1.2~1.8, 
with a mean value of ~1.5. Correlations have also been found 
between 'L2he.v, L2500 and aox, with the primary correlation 
being L2ksV oc L^soo- T^^^ slope index /3 was often found 
to deviate from unity for both optically selected (e.g. Strat- 
eva et al. 2005; Steffen et al. 2006; Just et al. 2007) and 
X-ray selected AGN samples (e.g. LussolO). However, La 
Franca, Franceschini & Cristiani (1995) re-analyzed Wilkes 
et al. (1994)'s sample by considering both variables and 
intrinsic scattering, and found that /? was consistent with 
unity. Green09 collected a large, well-defined sample of 2308 
SDSS/ChaMP QSOs in the redshift range 0.2~5.4, and con- 
cluded that P is not less than unity. 

The reality of a non-linear correlation in 'L2keV vs. 
L2500 remains an open question, but one possible explana- 
tion could be a selection effect in a flux limited sample for 
which dispersions in the optical and X-ray luminosity are not 
equal, or which change with cosmic time (Yuan, Siebert & 
Brinkmann 1998, hereafter: Yuan98; Tang07). However, the 
possibility of a truely intrinsic non-linear correlation cannot 
be ruled out. A non-linear L2kev vs. L2500 correlation implies 
that there is a dependence of aox on \j2kev and L2500 (e.g. 
Vignah, Brandt & Schneider 2003; Just et al. 2007; LussolO), 
but this is still a matter of debate (Yuan98; Tang07). 

To further test the basis of these correlations, we also 
calculated the values of L2keV , L2500 and aox from our best- 
fit model of the SEDs, and then preformed the same cross- 



correlation analysis. The limitation of our results arises from 
the lack of actual spectral coverage at 2500A for the 16 
sources without OM UVWl and UVM2 data. The lumi- 
nosity and redshift range of our sample is also relatively 
small. But on the merit side we have included two inputs of 
reddening/absorption to model both the Galactic and the 
AGN's intrinsic extinction during the broadband SED fit- 
ting, so our values of h2kev, L2500 and aox should be closer 
to those of the intrinsic source. The unobscured nature of 
our sample and the exclusion of warm absorber sources also 
helps reduce uncertainties in the corrections caused by red- 
dening/absorption. 

3.1.1 The L'ikeV vs. Z/2500 Correlation 

Figure [2] shows our L2feeV vs. L2500 correlation. The statis- 
tical methods used are the same as in LussolO, i.e. we use 
the full parametric estimate and maximized regression (EM) 
algorithm. We use this to derive two regression lines assum- 
ing, first h2kev, then L2500 to be the independent variable. 
Then the bisector of the two regression lines is calculated 
using the equations in Isobe et al. (1990). This method is 
more appropriate in cases where the cross-correlations are 
dominated by intrinsic scatter. The correlations found are 
as follows: 

(i) the EM regression line, when L2500 is assumed to be the 
independent variable: 

Log{L2kev) = (0.83±0.05)Lo5(L250o) + (1.59±1.55) (1) 

(ii) the EM regression line, when h2kev is assumed to be the 
independent variable: 

Log{L2kev) = (1.09±0.09)Lo(?(L25oo) - (6.17±2.66) (2) 

(iii) the bisector of the above two regression lines (the solid 
orange line in Figure [2]) 

Log{L2kev) = (0.95±0.06)Lo5(L250o) - (2.04±1.77) (3) 

The Spearman's rank test gives a rank coefficient of ps = 
0.87, and the probability of deviation from a random dis- 
tribution is ds = 1.2x10^^^, confirming a very high level 
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Figure 2. The cross-correlatin between L2500 arid ^2ke.v ■ The 
soUd orange line is the bisector regression line for our sample. 
The red triangle symbols represent NLSls; purple symbols show 
the radio loud AGN; the blue square symbol is Mrk 110. The filled 
purple star is the BAL-quasar PG 1004+130, and the open purple 
star is the position if its intrinsic X-ray flux was 0.73 dex higher 
(Miller et al. 2006). All triangle symbols represent Population A 
sources whose H/9 FWHM is less than 4000 km s~^. In the two 
histograms the green and red regions are for the Population A 
sources, and the red region indicates the 12 NLSls. The dashed 
green line is based on HasingerOS; the dash-dotted cyan line is 
based on Green09; the dotted pink line is based on LussolO. 



of significance. We superimpose PG 1004-1-130 (filled purple 
star) on Figure [2] showing that it lies far from the corre- 
lation due to its unusual X-ray weakness. It matches much 
better to the regression line if corrected in L2keV by 0.73 
dex (Miller et al. 2006, the empty purple star in Figure [Sjl. 

Our correlation between h2kev vs. L2500 is close to lin- 
ear, but previous studies have found a wide range of values 
as listed in Table 12.21 Figure [2] plots these results for com- 
parison. The correlation found by HasingerOS (green dashed 
line) which is based on an X-ray selected sample, is the 
most consistent with our sample, whereas the slope found by 
LussolO (pink dashed line) is significantly flatter. Our slope 
is also consistent with Green09 (cyan dashed line) in which 
a large sample of optically selected quasars is analyzed. We 
note that our sample only covers the low luminosity region 
of the sample in Green09. The NLSls are the least lumi- 
nous sources. The different value of the Y-axis intercept in 
Green09 may be due to their larger sample and larger lumi- 
nosity dispersion. 

There can be several reasons for the difference between 
our results and LussolO. Firstly there may be a selection 
effect of a flux limited sample if there are different amounts 
of dispersion in optical and X-ray luminosities (Yuan98; 
Tang07). More importantly, our sample is corrected for 
both Galactic and intrinsic reddening/absorption in the host 
Galaxy through the spectral fitting whereas that of LussolO 



is only corrected for Galactic absorption. We remove the 
intrinsic reddening correction and re-compute the EM re- 
gression, with results given in Table [22] under the row 'This 
Work-r^nt' The dust reddening and gas absorption column 
are related by, E(B-V) = 1.7x (NH/10^^)cm-2 (Bessell 1991), 
which means that L2500 is suppressed much more severely 
than \j2kev- Hence the removal of intrinsic reddening cor- 
rection decreases our correlation slope from 0.95 (±0.06) to 
0.91 (±0.05). Thus the intrinsic reddening can flatten the 
correlation, but it does not seem to be enough on its own to 
explain the difference with LussolO, unless their sample is 
strongly reddened in the optical. 



3.1.2 The aox vs. L2keV and L2500 Correlations 

To further investigate the correlation between the opti- 
cal/UV and X-ray continua, we adopt the same approach 
as in previous work to produce Oox vs. L2500 and aox 



correlations. If we assume L2 



« Lggoo, then 

aox oc L250Q and aox oc Lj^^^''^ are expected by deflni- 
tion. However, the Spearman's rank test does not imply very 
strong correlations: ps = 0.31,^3 = 0.03 for aox vs. L2500; 
Ps = — 0.13, ds = 0.35, for Oox vs. L2kev- The regression 
lines were derived but with large uncertainty. The results 
are presented in Figure [3] and listed in Table \2l2\ The solid 
orange line in each panel of Figure |3] is our EM regression 
line, compared with some previous work shown by dashed 
lines. Note that our results are limited to redshift 2 < 0.4, 
L2500 < 10+^' erg s-i Hz-i and LikeV < W+''^ '' erg s'^ 
Hz~^. Our results also suggest that the cross-correlations 
in aox vs. L2500 and L2keV are dominated by the AGN's 
intrinsic dispersion. 



3.2 Group 2: aox, ^2- 



and Xsdd 



As mentioned in the previous section aox is often used as 
a proxy for the broadband SED shape. Since h^oi is often 
dominated by the big blue bump (BBB) peaking in the un- 
observable EUV region (Walter & Fink 1993). K2-iofeeV, de- 
fined as \jboi /^2-iokeV , is also an indicator of the SED shape. 
XEdd is an important parameter which relates directly to the 
accretion processes close to the central SMBH. Therefore, 
correlations are to be expected between Xsdd, >^2-wkeV and 
aox • 



3.2.1 The aox vs. K,2-iokeV Correlation 

LussolO reported a tight second-order polynomial correla- 
tion for H2-wkeV vs. aox- We find a similar correlation for 
our sample. Spearman's rank test shows ps — 0.73 and ds 
= 2xl0~^. We also fitted a second-order polynomial to the 
correlation and obtained the following equation: 

Log{K2^io) = (5.7±3.4) - (8.8±5.2)ao, ± (4.3±2.0)aL (4) 

Figure [4] shows our best-fit polynomial (solid orange line) 
with ±l(j dispersion region (the shaded region). Note that 
our fit excludes BAL quasar PG 1004±130 (purple star in 
Figure |4)|. The best-fit polynomial from LussolO is plotted 
as the dashed green line, which is not as steep as ours. The 
reason is that our value of K2-iokeV (and Ltoi) is higher than 
found by LussolO, especially for the narrow line objects (the 
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Figure 4. re2-l0fceV vs. otox- Each symbol represents the same 
type of source as in Figure [2] The solid orange line is the best fit 
line found using a second order polynomial, and the shaded area 
is the itlfT zone. 



average K2-iofeev for our 12 NLSls is 86I45). LussolO con- 
structed their broadband SEDs by first assuming a power 
law extending from the optical to 500A, then connecting 
the continuum at 500A linearly to that at 1 keV, and finally 
by extrapolating from 1 keV towards higher energies, using 
an exponentially cut-off power law. This model substantially 
underestimates \jhoi for narrow line objects because such ob- 
jects often have strong soft-X-ray excesses which contain 
a large fraction of the Ljo; (Middleton et al. 2009; Jin et 
al. 2009; Paper-I). Our detailed broadband SED fitting has 



modeled this soft-excess feature by including a low temper- 
ature optically thick Comptonization component. We claim 
that this results in a more accurate estimate of h^oi (Paper- 
I). So certainly for the nearby Type 1 AGNs (redshift < 0.4), 
the K2-iofceV vs. aox correlation we find should be more re- 
liable. How the correlation behaves at high redshift requires 
further study, but LussolO has shown that such a second- 
order polynomial correlation still holds for Type 1 AGNs up 
to redshifts 2 = 4. 



3.2.2 The aox vs. Xsdd Correlation 

The existence of a correlation of Oox vs. XEdd remains un- 
clear. VF07 found no correlation between these quantities, 
and so they proposed that Oox did not provide useful infor- 
mation on the broadband SED shape. SOS confirmed VF07's 
result for their sample of 35 moderate to high luminosity 
radio-quiet AGN. On the contrary, LussolO did find a corre- 
lation between Uox and Xsdd, although with a large disper- 
sion. This was confirmed by GrupelO for their sample con- 
taining 92 soft X-ray selected AGNs, but their correlation 
was both fiatter and stronger than that of LussolO. We use 
our sample to investigate this situation, and our results are 
shown in Figure[5l The Spearman's rank test gives ps — 0.35 
and ds = lxlO~^, suggesting that a correlation does exist. 
We then applied the ordinary least squares (OLS) regres- 
sion, assuming Xsdd to be the independent variable, and 
found the following relation: 

aox = (0.079±0.038)Lo.g(Ai3dd) + (1.384±0.029) (5) 

In Figure[5]we show our results. Our correlation has less 
dispersion than found by VF07 and LussolO, but has larger 
dispersion than that from GrupelO. The exclusion of a corre- 
lation is at the '^2a significance level, which is less significant 
than in GrupelO. Our regression line slope is consistent with, 
but slightly fiatter, than that in LussolO (/3 = 0.133±0.023) 
and GrupelO (/? = 0.11±0.02). This is partly because our 
estimation of Ltoi is higher than in previous studies, due to 
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Figure 5. \Edd vs. ctox- Each symbol represents the same type 
of source as in Figure |2] The solid orange line is the bisector 
regression line determined by our sample. The cyan crosses are 
the binned data points of our sample. The dashed green line is 
from LussolO; the dash-dotted purple line is from GrupelO. 




Figure 6. A^^^ vs. K2— lOfceV- Each symbol represents the same 
type of source as in Figure [2] The solid orange line is the bisector 
regression line determined for our sample. The binned data points 
are from VF07 (pink) and VF09 (blue). The dashed green line is 
from LussolO. 



the inclusion of a soft X-ray excess in our model. Therefore 
our value of Xsdd is also higher for the NLSls and other 
relatively narrow line objects. Another reason could be a se- 
lection effect. We have shown in the previous section that 
aox has a luminosity dependence, thus at higher redshift 
we tend to detect more luminous sources with steeper aox- 
While both our sample and GrupelO's have low redshifts 
(us: z < 0.4; GrupelO: z < 0.3), LussolO's sample covers a 
much larger range in redshift (0.04 < z < 4.25), and indeed 
LussolO's sample contains many objects with aox > 1-5 and 
z > 0.4, which populate the empty region above aox ~ 1-5 in 
FigureO and create a larger dispersion. The large dispersion 
observed in our study and previous work suggest that one 
should be cautious about using the aox vs. Xsdd relation, 
because the Xsdd inferred by aox may contain considerable 
uncertainties. 



3.2.3 The K2-iokeV vs. Xsdd Correlation 

Wang, Watarai & Mineshige (2004) reported the correlation 
between K2-iofeeV and Xsdd, which was later confirmed by 
VF07,09. Most recently, LussolO also found this correlation 
for the 545 X-ray selected type 1 AGNs from the XMM- 
COSMOS survey. They suggested that the observed step 
change in this correlation does not result from the hboi de- 
pendence on both K2-iofceV and Xsdd- In our study we also 
find that Xsdd is correlated with K2-iofceV- A Spearman's 
rank test resulted in ps = 0.60, ds — 5xl0~^ for the whole 
sample, and ps = 0.60, ds = 5x10"" for the 12 NLSls. 
Figure [H] shows our results, together with the results from 
VF07,09 and LussolO. We performed an EM regression anal- 
ysis and derived the following equations: 



(i) An EM regression with Xsdd being the independent vari- 
able 

Log{K2-io) = (0.482±0.088)Loc/(AEdd) + (1.840±0.071) (6) 

(ii) An EM regression with H2~iokeV being the independent 
Log{n2-io) = (1.179±0.166)Lo3(Ai5dd) + (2.232±0.090) (7) 

(iii) bisector of the above two lines (solid orange line in Fig- 
ure [6]): 

Log{n2-io) = (0.773±0.069)Lo3(Ai5dd) + (2.004±0.049) (8) 

Our regression lines are highly consistent with the binned 
points from VF07,09 and also the regression line reported 
by LussolO. The two lowest data bins from VF07,09 seem to 
have a relatively high deviation from the correlation lines, 
which may imply a change in slope of the correlation for 
sources with Xsdd < 0.01. But we cannot test this possibility 
from our data due to the exclusion of sources with low Xsdd 
resulting from our sample selection (Paper-I). 

The results show that the K,2-iokeV vs. Xsdd correlation 
extends up to high K2-iofcey (~100) and super Eddington 
accretion rates (~10); such objects are mainly NLSls (red 
square symbols in Figure [6]) and some other relatively nar- 
row line sources (black square symbols). We also note that 
the dispersion in our regression line is smaller than that 
in VF07,09 and LussolO, in spite of the different methods 
used in deriving hboi and the different redshift ranges. This 
suggests that the dispersion observed in the correlation is in- 
trinsic. In Figure [61 we see that PG 1004-1-130 (filled purple 
star) deviates far from the regression line (also more than 3a 
from VF07,09's binned data points), confirming its anoma- 
lously weak L2_iofceV- Increasing its L2_iofeeV by 0.73 dex 
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Figure 7. K5100 vs. \Edd and K2-i0fceV- Each symbol represent the same type of source as in Figure [21 In the left panel, the solid 
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and \Edd=^- The symbols and lines in the right panel have the same meaning as those in the left panel. 



(open purple star) moves it much closer to the correlation 
line. 



3.3 Group 3: K5100, )^Eddr and K,2-iokeV 

The 5IOOA monochromatic continuum luminosity (L5100) 
is often used to estimate Ltoi, particularly for very large 
samples of AGN, when broadband SED modeling for every 
source is not practical. The conventional method is to use a 
constant scaling factor acsioo ~ 9 (Kaspi et al. 2000; Richards 
et al. 2006: 10.3±2.1), or «5ioo value that is anti-corrected 
with Lboi (Marconi et al. 2004, hereafter: Marconi04). VF07 
showed that for high XEdd sources such as many of the 
NLSls, there is a clear deviation from constant K5100 ~ 9. In 
addition, potential contamination from the host galaxy will 
introduce dispersion into the K5100 vs. XEdd correlation for 
low luminosity sources. However, this should not be a severe 
problem for our sample since in our sample host galaxy is not 
dominating (Paper-I). In our study a much stronger correla- 
tion was found in K5100 vs. XEdd as the Spearman's rank test 
gives ps=0.81 (ds=4x 10~^^). Motivated by the strong corre- 
lations between \Eddr and K2-iokeV, we also found a strong 
correlation between K5100 and K2-iofceV, with a Spearman's 
rank test of ps=0.64 (^3=9x10"^). 

The left panel of Figure[7]shows the correlation between 
K5100 and XEdd- The solid orange line is the OLS regression 
line, the two dashed orange lines show the ±lcr region, and 
the shaded region is the ±2(7 region. For a specific XEdd 
value, the la dispersion of K5100 is ~ 0.17 dex. The binned 
data points from VF07 are also shown in the plot for compar- 
ison. VF07's results are consistent with ours within ±2(t, but 
our correlation is much stronger. This may be attributed to 
the high spectral quality of our sample and the carefully de- 



rived Atsioo and Xsdd, based on our detailed broadband SED 
fitting. It also shows that the distribution of «5ioo peaks at 
10~20, with a la dispersion of 0.29 dex. For the NLSls, the 
mean K5100 increases to 20 (la — 0.23 dex). This means that 
using a K5100 = 9 (the horizontal purple line in the left panel 
Figure O would underestimate the intrinsic L(,oj and XEdd, 
especially for samples containing sources with high XEdd e.g. 
the NLSls. The OLS regression line that assumes Xsdd to be 
the independent variable can be expressed by the following 
equation: 

Loginsioo) = (0.467 ±0.045)Lo3(ABdd) + (1.430 ±0.027)(9) 

We superimpose the SDSS and FAST (Landt et al. 2011) 
data from Mrk 110 (filled and open green square, respec- 
tively) on Figure [T] This shows the large optical variability 
in the spectrum. The FAST data is much more consistent 
with the XMM-Newton OM and also matches very well with 
the regression line. This supports the reliability of the cor- 
relation. 

The strong correlation between K5100 and «:2-iofceV 
(shown in the right panel of Figure [7]) is an expected result, 
given that both K5100 and K2-iokeV strongly correlate with 
XEdd- Such a correlation is also reasonable considering the 
strong correlation between L2-iofceV and L5100 (Paper-II). 
The la dispersion of K2-iokeV is 0.40 dex, which is larger 
than that 0.29 dex for ksioo- The ±la, ±2a zones are also 
shown in the Figure. The bisector regression line can be ex- 
pressed by the following equation: 

Log{K5ioo) = (0.593±0.053)Loff(fi:2io) + (0.239 ±0.086) (10) 

Mrk 110 and PG 1004+130 are superimposed on the plot, 
and their corrected positions are much more consistent with 
the regression lines. 
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Marconi04 proposed that K5100 anti-correlated with 
hboi, but our study does not support such an anti- 
correlation, although our sample only occupies the Ljo; re- 
gion above 10^"'^ Lq in Fig. 3 left panel of Marconi04. 
A Spearman's rank test for our sample gives p3=0.12 
(ds=0.39), suggesting no correlation. The sources lying be- 
tween 10^^ ~10^^'^ Lq have a mean ^5100= 16, with a la 
dispersion of 0.29 dex. So we find that our K5100 values for 
these objects are much higher than reported in Elvis et al. 
(1994) and Marconi04. 

3.4 Group 4: r2-iofcey, ^Edd and i<i2-iokeV 

The strong correlation between 2-10 keV photon index 
{T2-wkev) and Xsdd has been studied in detail for the 
past ten years (e.g. Lu & Yu 1999; Wang, Watarai & Mi- 
neshige 2004; S06,08; ZhoulOa). It is proposed that increas- 
ing the mass accretion rate leads to enhanced emission from 
the accretion disc, resulting in more seed photons from 
the disc, which then increases the Compton cooling of the 
corona, and softens the Comptonized hard X-ray spectrum, 
i.e. the slope of V2-iokeV increases. It was also reported 
that both r2_iofceV and A^dd strongly correlate with the 
FWHM of (e.g. Brandt, Mathur & Elvis 1997; S06,08; 
GrupelO), therefore these three parameters all strongly cor- 
relate with each other. However, 806, 08 found that the cor- 
relation of V2-iokeV vs. FWHMh/3 is weakened by the in- 
clusion of highly luminous sources, but that the correlation 
of V2-iokeV vs. Xsdd Still exists. This implies that the cor- 
relation of T2-iokeV vs. Xsdd is morc fundamental. We also 
mentioned in Section[3]2]that the strong correlation between 
Xsdd and K2-iofceV is confirmed, thus a strong correlation 
between V2-iokeV and H2-wkeV is expected. Indeed, such 
a correlation has been reported recently by Zhou & Zhao 
(2010), hereafter: ZhoulOb. In this section we carry out a 
similar cross-correlation study, to test the robustness of pre- 
vious claims. 



3.4-1 The Correlations and Regression Lines 

The two panels in Figure [S] show our correlations between 
^2-wke.v, Xsdd and K2-iofceV. Table [3] Summarizes the nu- 
merical results. Significant correlations are confirmed based 
on the Spearman's rank test: p3=0.40 (^3=4x10"^) for 
r2-iofcev vs. Xsdd, and ps=0.73 (ds=4xl0"®) for V2-iokev 

vs. K.2-U)ke.V ■ 

Following S08's approach, we applied the minimiza- 
tion method for T2-\okeV vs. Xsdd correlation, assuming 
r2-iofceV = PLog{XEdd) + A typical error of 10% was 
assumed for Xsdd- The small error in I'2-iokeV for Mrk 
110 (the blue square symbol) caused the slope jS to be 
0.018±0.019, which is clearly not the best-fit line for the 
whole sample. We therefore excluded Mrk 110 and so found 
a more reasonable slope of 0.189±0.026, but this is still 
~5a away from 0.31±0.01 reported by SOS using the same 
method. It implies that the minimization technique may 
not be an appropriate method for quantifying this correla- 
tion, because it can be strongly biased by sources with small 
error in the r2_iofeeV measurement (if the 2-10 keV spec- 
trum has high S/N). The x^/*^ — 6.5 in our fitting means 
that this correlation contains a big intrinsic dispersion along 



Table 3. The line coefficients found using different regression 
methods for the correlations of r2-iofceV vs. Ag^jd a-nd i^2—l0keV- 







Log{XEdd)=l3rH 


iog(K2io)=/3r-f,e 


OLS(Y|X) 




0.918±0.269 


1.115±0.172 




? 


-2.325±0.527 


-0.573±0.329 


OLS(X|Y) 


13 


4.730±1.217 


2.209±0.283 




? 


-9.650±2.304 


-2.675±0.525 


Bisector 




1.753±0.239 


1.533±0.153 




? 


-3.931±0.471 


-1.376±0.288 


min 




4.416±0.510 


2.620±0.184 




? 


-8.761±0.946 


-3.447±0.343 






1.274±0.283 


1.529±0.183 






-3.003±0.544 


-1.372±0.350 



with the observational dispersion, thus the method of assum- 
ing /'^ ~1 by taking intrinsic dispersion into account is 
more appropriate. This method gives slopes of 0.202±0.061 
and 0.226±0.026 before and after excluding Mrk 110, so 
the results are less sensitive to individual sources of much 
smaller error bars. The intrinsic dispersion is 0.18, which 
is 86% of the total dispersion, and is also consistent with 
-lokev ^0.1xF2-iokeV reported in SOS. The bisector 
regression line is derived. The result is plotted as a solid or- 
ange line in the left panel of Figure[Sl and it can be expressed 
by the following equation: 

Log{XEdd) = (1.773 ± 0.23S)r2-iofcev - (3.9S3 ± 0.469)(11) 

The slope is ~2(t steeper than the bisector slope of 0.9±0.3 
reported by SOS (dashed green line in the left panel of Fig- 
ure|8]). This discrepancy may be due to the different methods 
used to estimate the bolometric luminosity. We will discuss 
this point in Section [3.4.21 

Similar analytical methods were applied to the rela- 
tion of T2-iokeV vs. K2-iokeV ■ ZhoulOb reported a slope 
of 2.52±0.08, using standard x^ minimization and assum- 
ing K2-iokeV = l3Log{r2-iokev) + Thls is consistent with 
our value of /3 = 2.620±0.1S4 with a xV'^ = 2. 78. Con- 
sidering the intrinsic scatter, ZhoulOb reported a slope of 
1.12±0.30 by adding 0.32 dex of intrinsic dispersion to re- 
duce x^/i^ to unity (the dashed green line in the right panel 
of Figure [8]) . Applying the same method to our sample re- 
sulted in a slope of 1.529±0.183, which is steeper than found 
by ZhoulOb. The intrinsic dispersion found by us is 80% of 
the total dispersion. The bisector regression method gives a 
slope of 1.533±0.153 (the solid orange line in the right panel 
of Figure [S]) and can be expressed by the following equation: 

Log{K2-iokev) = (1.533 ± 0.153)r - (1.376 ± 0.288) (12) 

This is also consistent with the results found by assuming 

xV^ -1. 

3.4-2 Advances from Our Correlation Analysis 

Compared with the results found by S08 and ZhoulOb, 
our study of Fa-io^eV vs. Xsdd and K2-iokeV provides the 
following advances. 

(1) We have confirmed these correlations based on 
sample of twice the size of those in S08 and ZhoulOb, 
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Figure 8. lofcev vs. XEdd and K,2~iokeV ■ Each symbol represents the same type of source as in Figure[2] In the left panel, the solid 
orange line is the bisector regression line. The dashed green line is that reported by ZhoulOb. The vertical and horizontal purple lines are 
for r2_iofceV=2 and \Edd=^- The symbols and lines in the right panel have the same meaning as those in the left panel. The vertical 
purple line is for ft2— iOfceV=100. 



including more sources with high values of Xsdd, ^2-iokeV 
and which significantly extend the previous 

correlations (see Figure [8|. The regression hne fits are 
better constrained and cover wider parameter space. The 
difference between our regression lines and those of previous 
studies may be partially due to the fact that we have more 
sources of most extreme Xsdd- 

(2) Our sample has been carefully screened to exclude 
sources with a strong warm absorber. These sources may 
have r2-ios;e\^ and higher K2-iokeV artificially lower than 
the intrinsic values. Our sample quality is essential to 
reduce the dispersion and so reveal intrinsic correlations. 

(3) Our estimates of Ltoi were derived from the broadband 
SED fitting, which was based on high quality spectra 
and a new multi-component model. We claim this to 
be more reliable than the procedure used in previous 
studies. A conventional method is to apply a multiplication 
factor to L5100 to estimate Ltoj. However, we showed in 
Section [3.31 that K5100 is well correlated with Xsdd, rather 
than being constant or dependent on L(,oj, consequently 
the conventional scaling from L5100 to Ltoi is likely to 
result in poor accuracy in some cases. The Lboi used in 
ZhoulOb does come from VF09's broadband SED model 
for the reverberation mapped sample, but it does not take 
account of the 'soft X-ray excess' component or where the 
disc peaks in the EUV. Therefore the L^oi we calculate will 
be larger than previous works, especially for those sources 
with a strong 'soft excess', Our Xsdd and K2-iokeV will also 
be higher, which could account for the differences in slope 
between our regression lines and those reported in previous 
work. 



3.5 Group 5: r2-iofceV, H/3 FWHM and Mbh 

3.5.1 The V2-iokeV vs. H(3 FWHM Correlation 

The correlations between the soft/hard X-ray slopes and 
the Balmer line velocity width have been the subject of 
many previous studies. Puchnarewicz et al. (1992) studied 
53 AGNs with ultra-soft X-ray excess and noticed that these 
ultra-soft AGNs tend to have narrower optical permitted 
hues. Laor et al. (1994) studied 23 ROSAT selected bright 
quasars, and found an anti-correlation between the 0.2-2 keV 
slope {ax) and the H/3 line width. Later Boiler, Brandt & 
Fink (1996) showed that NLSls tend to have softer X-ray 
spectra. Brandt, Mathur & Elvis (1997) extended this anti- 
correlation to include the 2-10 keV slope, by showing that 
NLSls also have steeper hard X-ray continua than BLSls, 
a result which was confirmed and extended by other stud- 
ies (e.g. Grupe et al. 1999; Leighly 1999; Piconcelli et al. 
2005; Brocksopp et al. 2006; SOS; ZhoulOa). However, there 
is large scatter within this correlation, and the trend seems 
to invert below ~1000 km s"^ for NLSls (Zhou et al. 2006). 
The observed large scatter is to be expected since we know 
that it is not a single variable that determines the spec- 
tral slope in either the soft or the hard X-ray region. In the 
soft X-ray region, the extinction, a soft-excess component 
and a warm absorber will all infiuence the spectral shape, 
which would require very detailed modeling. The situation 
for the 2-10 keV region is somewhat less complicated since 
often a single power law dominates (e.g. Middleton, Done 
& Gierliiiski 2007) , but a warm absorber and refiection may 
still modify the hard X-ray spectral shape. In summary, the 
true correlation can only be found when the intrinsic X-ray 
continuum is used. 

Our sample selection has ensured that every object in 
the sample has high quality 2-10 keV spectra, without sig- 
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Figure 9. lOfeeV vs. FWHM of H/3 and M^^f . In the upper left panel, a broken line is fitted to the sample using the minimum 
method, with the break point chosen to be FWHM = 4000 km s~^ (the orange square point). S06,08 proposed the linear correlation 
between r2_iofeeV a-iid Log(H/3 FWHM) was not followed by their 10 extremely high luminosity sources, so wc plot their sample as blue 
diamond symbols for comparison. In the right panel binned points are plotted with 1 standard error of r2_iofeeV in order to show the 
break points more clearly. The two red points only include the NLSls, the two dark points are the broadest BLSls. The cyan point is 
the binned point for the whole sample of S06,08. The blue point is the binned point for S06,08's sample but excluding LBQS 0109+0213 
whose r2_iofceV is anomalously low. IE 1556+27.4 shown by the red circle is another source with r2_iofeeV <1.5. In the second row, 
the symbols all have the sample meaning as in the first row, the break point is chosen to be Log(Ma£f) = 8.0 (the orange square point). 
We plot the linear regression line as the dashed orange line. 



nificant cold gas absorption or a warm absorber (Paper- 
I). We confirm that there is an anti-correlation between 
r2-iofeeV and H/3 FWHM, see Figure [9] The Spearman's 
rank test gives pg =-0.72 (ds=4.9xl0~^). The best-fit lines 
from ZhoulOa are also plotted in Figure [5] as the cyan lines. 
Compared with their linear correlation in using FWHMh^i, 
we find that the linear correlation using Log(FWHMH^() is 
better. Previous work also noted that the correlation may 
change form at ~4QQQ km s~^ (Sulentic et al. 2008). There- 
fore we fit a broken power law to the data points in the 



left panel by assuming a break point at FWHMff^ = 4000 
km s~^. Intrinsic dispersion will dominate this correlation as 
shown in previous studies (e.g. Grupe et al. 1999; GrupelO; 
ZhoulOa), so we performed this fitting without considering 
the error in V2-\okeV associated with each point, otherwise 
the the fitting would be biased by the few points with best 
constrained r2-iofeeV rather than revealing the distribution 
of the whole sample. The best-fit parameters are shown be- 
low: 

(i) when FWHMh/3 4000 km s~\ 
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Figure 10. Mbh vs. H/3 FWHM. Note that the black hole mass 
M-BH was constrained by H/3 FWHM and obtained from the 
broadband SED fitting (i.e. the Mbh FIT in Paper-I), rather 
than directly calculated from the H/3 FWHM. The symbols used 
represent the different types of source as in Figure [2] The solid 
orange line is the OLS regression line, assuming H/3 FWHM to 
be the independent variable. The shaded region is the itlo" re- 
gion of the regression line. The cyan triangle shows the position 
of (FWHMjf/3,6reafe, Log(MB if ,6r-eafe )) m Figure[9l 



r = (-0.86 ± Om)Log{FWHMHii) + (4.89 ± 0.03) (13) 
(ii) when FWHMh,? > 4000 km s"\ 

r = (-0.11 ± Om)Log{FWHMHii) + (2.20 ± 0.04) (14) 

bources with FWHMif,3 > 4000 km s"^ have an average 
r2-iofceV = 1.78±0. 12. The only source included in the cor- 
relation whose r2-iofceV < 1-5, is IE 1556-1-27.4, (the data 
for PG 1004-f 130 is superimposed but not used for the re- 
gression). A closer examination of the spectrum of this AGN 
shows that it probably has a strong reflection component 
modifying its intrinsic hard X-ray power law slope (Paper- 
I). All other objects have values consistent with r2_iofeeV > 
1.5. The differences in the results of ZhoulOa and our work 
are not only because we performed our correlation fitting us- 
ing Log(FWHM£f^), but also because their sample did not 
exclude BAL quasars and warm absorbers, whose low val- 
ues of r2-iofceV are probably not intrinsic. This will bias the 
correlation and increase the scatter. 

Our sample includes six objects with FWHM///3 > 10000 
km s~^. These are the sources have average < r2~iofceV >= 
1.76±0.14 independent of the FWYiMap- This is slightly 
lower but still consistent with the r2-iofceV = 1.97 ± 0.31 
index found by S06,08, who included more high redshift, 
high luminosity sources, with FWHMif/j > 10000 km s"^ . 



3.5.2 The V2-iokeV vs. Mbh Correlation 

The H/3 FWHM is frequently used to estimate the Mbh, 
using the relation M.bh oc FWHM|j^^ (Wandel, Peterson 
& Malkan 1999; Woo & Urry 2002). The correlation of 
r2-iofceV VS. FWHM^f/j implies a similar correlation in 
r2-iofceV vs. Mbh- This is confirmed in our study as shown 
in the second row of Figure [9] The Spearman's rank test 
gives ps=-0.3 (ds=3xl0"^). We also plot FWHM///3 vs. 
Mbh in Figure [TOl This is an independent plot as our Mbh 
are derived from the SED continuum fits rather than di- 
rectly measured from FWHMff/j. The OLS regression gives 
Mbh tx EWnM]^^**, with a Icr = 0.4 dex. The cyan trian- 
gle symbol in the figure shows the position of the break in 
FWHMff/j vs. Mbh. The correlation in Figure [TOl suggests 
FWHMif/3 = 4000 km s~^ corresponds to a black hole mass 
of 10* Mq. Indeed, the correlation in the second row of Fig- 
ure[5]shows a change in slope at Log(Msi/) ~ 8.0. A broken 
power law fitting, assuming the break point at log(Ms£r) = 
8.0, can be expressed as follows: 

(i) when Log(MBij) ^ 8.0, 

r2-iofcei' = (-0.372±0.005)Log(MsH) + (4.802±0.037)(15) 

(ii) when Log(MsH) > 8.0, 

r2-iofcev = (0.056 ± Q.m7)Log{MBH) + (1.380 ± 0.052)(16) 

Sources with Log(MsH) > 8.0 have an average 
V2-10keV = 1.86±0.21. 



3.6 Other Strong Correlations 

The diverse correlations found in small parameter groups 
imply more correlations among all these parameters. We 
show some significant correlations in Figure 1121 These in- 
clude decreasing XEdd, K2-iokeV , K5100 with increasing H/3 
FWHM and Mbh. Binned data points are shown as cyan 
crosses. The dashed orange line in each panel is the bisector 
regression line. 



3.6.1 The K2-iokeV vs. Z/2-iofceV correlation 

Marconi04 and Hopkins07 reported a strong positive cor- 
relation in K2-iofceV vs. L2-iofceV, bascd on a quasar SED 
template and the aox vs. L2500 correlation reported by Stef- 
fen et al. (2006). VF07 tested the same correlation in their 
low redshift AGN sample but found no correlation (see Fig- 
ure [11]). We tested this correlation in our sample and con- 
firmed VF07's result (see Figure llip . A highly dispersed 
anti-correlation of K,2-iokeV vs. lj2-iokeV is found in our 
study. 

For the well sampled high X-ray luminosity region 
(L2_iofceV >10*^ erg s~^), both VF07 and our samples show 
very similar distribution and dispersion. NLSls in both sam- 
ples mainly populate the regions of K2-iofeeV -c^lOO and 10*^ 
erg s~^ <L2-iofceV <10*'^ erg s~^, and deviate from the cor- 
relation suggested by Marconi04 and Hopkins07. The main 
difference from VF07 is that our sample have more sources 
within the range of 3x10''^ erg s~^ <L2-iofceV <10''^ erg 
s~^, while VF07 sample contains three extraordinary weak 



X-ray sources whose L2_iofeeV < 2x10*^ erg s-\ Thus 
we think those intrinsically X-ray weak (L2-iofeeV ~10*^ 
ergs~^) sources may populate the low L2-iofeeV, small 
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K2-iokeV region, creating a totally dispersed distribution 
in the K2-iokeV vs. L2_iofceV plot. The correlations from 
Marconi04 and Hopkins07 may have underestimated the 
uncertainties in using aox vs. L2500 correlation (see discus- 
sion in Section rs.ip and the universal quasar SED template 
(e.g. Elvis et al. 1994; VF07; Paper-I). However, we cannot 
rule out the possibility that the behaviors of nearby Seyfert 
AGNs are different from quasars at higher redshift. 



3.6.2 otuv and ax related correlations 

While r2_iofceV is the 2-10 keV photon-index measured di- 
rectly from the X-ray data, we can also measure the opti- 
cal/UV slope and soft X-ray slope from our reconstructed 
SED. Assuming Vi, oc u~" , we define auv as the spectral 
slope between 1700-6500 A (the Swift UVOT wavelength 
coverage), and ax as the soft X-ray slope between 0.2-2 keV, 
so as to be comparable with the results in GrupelO. Note 
that we do not have complete OM data for every source, and 
we only use X-ray data from XMM-Newton above 0.3 keV, 
so our auv and ax values are model dependent. The soft 
X-ray and optical/UV regions cannot be simply fitted by a 
single power law (see the SEDs in Paper-I), and thus auv 
and ax are just rough estimates of the spectral shape. We 
measure auv and ax from the best-fit SED in Paper-I cor- 
rected for Galactic and intrinsic reddening/absorption, and 
list the values in Table [B2l 

It was found previously that AGNs with bluer opti- 
cal/UV spectra tend to have softer X-ray spectra (Wal- 
ter & Fink 1993; Grupe et al. 1998; GrupelO). We found 
ps = —0.65 {ds — 4x10"^) for auv vs. r2-iofecV, and 
ps — —0.41 {ds — 3x10"'^) for auv vs. ax, so confirm 



the results from other studies. We also tested the results 
of auv and ax without correction for intrinsic redden- 
ing/absorption. The Spearman's rank coefficients remain 
similar which is because of the unobscured nature of our 
sample. Another correlation found in previous work is that 
AGNs with steeper X-ray spectra tend to be weak at hard 
X-ray energies (Atlee & Mathur 2009; GrupelO; but see also 
Young, Elvis & Risaliti 2009). This can be directly confirmed 
in our study by the anti-correlation found between ax and 
^2-iokeV {ps ~ —0.60, ds = 6xl0~®), and the correlation 
between ax and aox {ps = 0.65, ds = 5xlO~^, Figure [13]). 
However, we only find a marginal anti-correlation between 
auv and aox (see Table [4|, which is much less significant 
than found by GrupelO. This seems to be due to the fact that 
our sample has few sources with aox > 1.6 and auv < 0. 

Similar to the T2-iokeV vs. XEdd correlation discussed in 
Section r3.4l ax and auv were also found to (anti-)correlate 
with XEdd (Grupe 2004; SOS; GrupelO). This is confirmed 
in our study as we find ps = 0.49 {ds = 4xl0~*) for ax 
vs. Xsdd which is stronger than the correlation between 
r2-iofcev and Xsdd, and ps = -0.55 {ds = 4x10"*) for auv 
vs. XEdd. In addition, we also find strong (anti-)correlations 
such as ax (or auv) vs. K2-iofeeV, ax (or auv) vs. K5100A, 
ax (or auv) vs. FWHM^f^ and ax (or auv) vs. Mbh (see 
Figure [T3l Table g]). 



4 A SYSTEMATIC CORRELATION STUDY 
ON THE KEY PARAMETERS 

To summarize the various correlations discussed in the pre- 
vious section, we performed a systematic correlation study 
of the following key parameters: r2_iofceV, «2-iofceV, K5100A, 
Xsdd, FWHMh,9, Mbh, aox, Lboi, ^'2-iokeV, auv and qjc. 

First, a correlation matrix was constructed as shown in 
Table |4l The Spearman's rank coefficient and probability of 
a null hypothesis are given. The table shows that there are 
some sub-groups of parameters which are strongly coupled 
with each other. For example, H2-iokeV , i^siooA and Xsdd 
are coupled; r2-iokev, H/3 FWHM, ax and K2-iokev are 
also coupled with each other. The strong correlation between 
H/3 FWHM and K2-iokeV can be expressed by the following 
equation derived from a bisector regression analysis: 



Log{n2- 



(-1.22±0.12)Log(Fl¥J/M)4-(5.88±0.45)(17) 



There also appears to be a sub-group consisting of Mbh, 
Lboi and L2-iofceV, and a sub-group consisting of Mbh, H/3 
FWHM and L2_iofceV. However, correlations within these 
sub-groups are probably a result of selection effects arising 
from our sample selection criteria. For example, he inclusion 
of extremely weak L2_iofceV sources may weaken or elimi- 
nate the correlations between L2-iofceV and other parame- 
ters. 

The observed properties of AGN should be ultimately 
driven by the black hole mass, mass accretion rate, black 
hole spin and orientation angle. We have assumed the sim- 
plest Schwarzschild black hole in our model and so its spin is 
not considered. Uncertainties introduced by orientation an- 
gle should also be small since our sample only contains un- 
obscured Type 1 AGNs. Therefore, the remaining intrinsic 
parameters are just the black hole mass and mass accretion 
rate (or equivalently, Eddington ratio) 
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Figure 12. Examples of some good correlations not reported previously. In each panel the various symbols represent the same types of 
source as in Figure [2] The cyan symbols are the binned data points over the X-axis with a 1 standard error on the Y-axis. 



A 

A 

• ★ 


A 

. ^'^ 

'ikti A 
*. ★ 


A 

★ • 


A 
A 

. ★ 


A 

^ * 


• 

A 

• 

A^r 

A 


• 

A 

A 


4-' 

• 

* A 

'^^V»A ■ . 
A 


• 

A 


• 

r*i A ^ ★ 

A 



0.01 0.10 1.00 10.00 10 



100 1000 10 100 1000 100001.0 1.2 1.4 1.6 1.8 2.0 

-lokev fsiooA H/? FWHM 



Figure 13. Examples of some good correlations not reported previously. In each panel the various symbols represent the same types of 
source as in Figure (2] The cyan symbols are the binned data points over the X-axis with a 1 standard error on the Y-axis. 



We can examine the correlations further by performing 
a principal component analysis (PCA) on the correlation 
matrix formed by Pearson's correlation coefficient (Pearson 
1901; Boroson & Green 1992; Francis & Wills 1999). First, 
we include all the 11 parameters and so the dimension of the 
correlation matrix is 11. Therefore the outcome of the PCA 
must contain 11 normalized eigenvectors (principal compo- 



nents: PCs), each associated with a positive eigenvalue. Each 
PC is a linear combination of the 11 parameters, and is or- 
thogonal to all the other PCs. The sum of the 11 eigenvalues 
equals 11. A higher eigenvalue would suggest a larger frac- 
tion of correlations contained in the direction of the corre- 
sponding eigenvector. 

The EIGENQL program in IDL (Interactive Data Lan- 
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Table 4. The cross-correlation matrix of the 11 key parameters, ps is the Spearman's rank coefficient, d], is the logarithm of the 
significance level of being random distribution. 
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log 
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log 
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1 


0.73 


0.32 


0.40 


-0.72 


-0.33 


0.39 


0.05 


-0.38 


0.63 


-0.41 




dl 


-oo 


-8. 


-2. 


-2. 


-8. 


-2. 


-2. 


-0. 


-2. 


-6. 


-3. 


K-2-lOkeV 


ps 


0.73 


1 


0.64 


0.60 


-0.81 


-0.45 


0.74 


0.12 


-0.49 


0.89 


-0.59 




dl 
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-0. 
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0.32 
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0.64 
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0.78 
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0.78 


-0.17 
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ps 
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guage) v6.2 was used to perform the PCA. We found that 
the first three eigenvectors contain 87% of the total corre- 
lations in the matrix, i.e. 53% in eigenvector 1 (PCI), 25% 
in eigenvector 2 (PC2) and 9% in eigenvector 3 (PC3). To 
determine the actual contributors of these eigenvectors, we 
cross-correlated them with the 11 key parameters. Table [5] 
lists the Spearman's rank coefficients. It is clear that PCI 
correlates/anti-correlates very well with most key parame- 
ters, except for Lboi- The highest correlation strength is for 
K2-iofeei/, ctx and H/3 FWHM. These in turn are probably 
proxies for the physical variables of Mbh and \Edd- PC2 is 
dominated by Lboi which confirms that Lboi is a relatively 
independent variable. We have assumed that Ltoi ~ fiMc^, 
where ^ is the standard accretion efficiency of 0.057 (see 
Paper-I), so PC2 is in effect dominated by the mass accre- 
tion rate M. The contribution from the rest eigenvectors to 
the set of correlation is small compared to PCI and PC2, 
and are therefore not important. It was also reported by 
Boroson (2002) that the PCI from the correlation matrix 
of optical emission line parameters is driven predominantly 
by the Eddington ratio, while the PC2 from the same ma- 
trix is dominated by the luminosity. Therefore, the PCI and 
PC2 from our correlation matrix of the 11 parameters have 
a similar basis to the two eigenvectors reported by Boroson 
(2002). 

The inclusion of \Edd, M_bh and Lboi may bias the out- 
come of PCA, so we performed an independent PCA by 



excluding XEdd, Mbh and Lboi from the correlation matrix 
(Grupe 2004; Grupe 2011), and so the remaining matrix only 
has a dimension of 8. The resultant eigenvectors were again 
correlated with each of the 11 parameters. The results are 
listed in Table [5] The two new principal eigenvectors (PCl- 
excl and PC-excl in Table [5]) are quite similar as PCI and 
PC2 in terms of correlation outcomes with the 11 parame- 
ters. This confirms our previous PCA finding that the three 
physical parameters Xsdd, Mbh and Lboi drive the majority 
of the correlations. 

5 THE MEAN SEDS 

5.1 Diversity of the Mean SEDs 

In Paper-I we derived the mean SED for the 12 NLSls in 
our sample, and compared this with the mean SEDs of an- 
other two groups of BLSls. We showed that the SED shape 
changes dramatically as the H/3 FWHM increases. Since 
quasar SED are not uniform, modelling their spectra and 
evolution using a single template such as that of Elvis et 
al. (1994), will mask out a large dispersion in their intrinsic 
properties. In our study we aim to minimize this dispersion 
by grouping the SED based on each of the 11 key parameters 
discussed in the previous section. However, since auv and 
ax are rough estimates of the spectral shape in optical/UV 
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Table 5. The cross-correlation of the eigenvectors with the 11 key 
parameters. The coefficients given are from the Spearman's rank 
test. PCI and PC2 are the principal components from the PCA 
on all the 11 parameters, while PCl-excl and PC2-excl are the 
principal components after excluding \Edd) ^BH and ^bol from 
the PCA. The row 'Property' shows the percentage of correlations 
contained in the direction of that eigenvector. 





PCI 


PC2 


PCl-excl 


PC2-excl 


Property 


53% 


25% 


62% 


14% 


r2-10fceV 


0.71 


0.03 


0.75 


-0.14 




0.93 


0.14 
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0.06 
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0.74 


0.26 


0.67 


0.34 


^Edd 


0.60 


0.67 


0.51 


0.60 


FWHM 


-0.89 


0.17 


-0.90 


0.21 


Mbh 
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0.48 


-0.65 


0.34 


aox 


0.59 


0.36 


0.63 


0.35 




-0.06 


0.99 


-0.10 


0.81 


L2-10fceV 


-0.60 


0.76 


-0.65 


0.67 


ax 


0.92 


-0.03 


0.93 


-0.01 


O-VV 


-0.74 


-0.04 


-0.70 


-0.17 



and soft X-ray regions, we excluded them from the mean 
SED calculation. 

There are 51 AGNs in our sample. For each of the 9 re- 
maining parameters, we sorted the sources according to the 
parameter value, and then classified the sample into three 
subsets so that each subset includes 17 sources. Then the 
BAL quasar PG 1004-1-130 was excluded from its subset. 
The SEDs constructed using Model-B in this paper were first 
corrected for redshift, and then divided into 450 energy bins 
between 1 eV and 100 keV in the logarithm. Within each 
subset we calculated the average luminosity in every energy 
bin in logarithm, together with the Icr deviation. Then we 
derived the mean SED for each group together with the la 
dispersion. The same calculation was repeated for each of 
the 9 parameters, so that there are three mean SEDs for 
each parameter. No special note was made for the NLSls 
because their defining parameter, H/3 FWHM, is one of the 
9 parameters. 

Figure[T3]shows our results. Each row displays the three 
mean SEDs divided according to the parameter shown in 
the panel title. The SEDs have all been renormalized to 
the mean luminosity at 2500A of each subset. To highlight 
the differences among these SEDs, we mark the locations of 
2500 A and 2 keV by the vertical orange lines. The relative 
height of these two lines directly reflects the value of aox, and 
the height of the bar at 2 keV shows the dispersion in 
since the co-added SED's are all normalized at 2500 A. We 
also mark the energy peak for each mean SED by the vertical 
thick purple line. The mean values of other parameters are 
given in each plot for comparison. 

We find that the mean SED changes in a similar way 
with all the parameters except 'Lhoi- The energy at which 
the disc emission peaks decreases along with the ratio of lu- 
minosity in the disc compared to the Comptonised compo- 
nents, and the 2-10 keV spectral slope. Our stringent sample 
selection means that these spectral differences are intrinsic 
rather than due to absorption/extinction. 

If the SED changes are determined solely by one of the 
parameters considered here, then binning based on that pa- 
rameter should result in the minimum dispersion within each 



individual binned SED, and maximal difference between the 
three SEDs spanning the range in that parameter. However, 
the SED changes should ultimately be physically dependent 
on changes in Msij and \Edd- Since there are two depen- 
dent variables no single parameter will completely deter- 
mine the behaviour. Hence the dispersion within each of 
the three binned SEDs is minimized (and the difference be- 
tween them is maximized) for composite parameters which 
depend on both \Edd and Mbh such as «;2-io/seV, K5100 and 
H/3 FWHM rather than the fundamental physical parame- 
ter \Edd- Future work with larger samples can improve on 
this study by selecting a subsample of AGN with different 
Afidd but similar Msh. Such mass selected samples would 
give the best comparison to the SED changes in Black Hole 
Binaries (BHB), which all have the same mass to within a 
factor of ~ 2. 

5.2 Discussion of the Mean SEDs 

Since our SED fitting is based on a physical model, we can 
'correct' for the mass dependence of the SED shape to get 
an estimate for the SED differences in AGN as a function 
solely of \Edd- This is shown in Done et al. (2011) for a 
Mbh ~ 10* Mq and forms the basis of a direct comparison 
with the BHB spectral states seen for these Mbh = IOMq 
systems as \Edd changes. This has many superficial simi- 
larities to the dramatic state change seen in BHB as their 
luminosity increases. The SED changes from a 'low/hard 
state' being dominated by Comptonisation, with a hard X- 
ray spectral index r2-io*:eV < 2, and the disc component 
peaking at rather low temperature, to a 'high/soft state' 
where the disc dominates the luminosity and the X-ray spec- 
tral index is softer, T2~iokeV ~ 2 — 2.2 (see e.g. the review 
by Done, Gierlihski & Kubota 2007). However, this occurs 
at Xsdd ~ 0.02 for moderate changes in mass accretion 
rate (Maccarone 2003), an order of magnitude lower than 
the spectral differences seen here in the AGN (see also the 
discussion in Done et al. 2011). This could indicate some 
subtle differences in the transition due to the very different 
masses, but in BHB this transition is also associated with 
the radio jet switching off (Fender, Belloni & Gallo 2004). 
If the Compton dominated states in AGN correspond to the 
low/hard state in BHB then we would expect them to be 
radio loud. However, AGNs are radio-quiet by a factor of 
10:1. 

This makes unlikely an identification of the two lower 
Asdd AGN templates as analogs to the low/hard (or inter- 
mediate state) seen in BHB. Instead, there is another state 
in BHB called the 'very high' or 'steep power law state', 
where the disc also peaks at a lower temperature than ex- 
pected, and where the Comptonised component contains a 
large fraction of the total luminosity (see e.g. the review by 
Done, Gierlihski & Kubota 2007). However, this state has 
r2-iofceV > 2.5 i.e. the hard X-ray spectra are steep. Yet 
these AGN have r2-iofeey < 2, as well as a soft X-ray ex- 
cess component. Therefore, to match the AGN with the very 
high state would require that refiection and/or complex ab- 
sorption modify the spectrum, producing an apparent soft 
X-ray excess and a hard power law from an intrinsically 
steep spectrum. However, the time variability properties of 
individual objects make it clear that these high mass, rela- 
tively low Eddington ratio objects do indeed have two sepa- 
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Figure 14. The AGN mean SEDs based on different values of the 9 key parameters from Model-B fitting (i.e. including the effect of a 
colour temperature correction). For each parameter, the 51 sources are sorted according to the parameter value, and then are divided 
into three equal subsets so that each contains 17 sources. PG 1004+130 is excluded from its subset. Finally, a mean SED is constructed 
for each of the three subsets after renormalizing each individual SED to the mean luminosity at 2500A of that subset. The three panels 
(A, B, C) in each row show the mean SEDs for the subsets classified by the parameter shown in the panel title. In each panel the solid 
curve is the mean SED, while the shaded coloured region is the itlcr deviation. The 2500 A and 2 keV positions are marked by the 
vertical solid orange lines, whose relative height indicates the value of ciox- The peak position of the SED is marked by the vertical solid 
purple line. The average values of some other parameters in that subset are also shown in the panel. Each mean SED has been rescaled 
by the same arbitrary constant on the Y-axis which is 1x10"''®. Note that the energy ranges 1.4 eV < E < 6 eV and 0.3 keV < E < 10 
keV, are covered by SDSS, OM and EPIC data respectively, while the SED in the rest energy bands is determined from an extrapolation 
of the best-fit model. 
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rate components. The intrinsically hard power law which is 
more variable on short timescales, and a soft X-ray excess 
which is relatively constant on these timescales but more 
variable over longer times (e.g. Mkn509: Mehdipour et al. 
2011; Noda et al. 2012). 

Therefore, the two lowest Xsdd spectra shown in Fig- 
ure [14] do not look similar to any state observed in BHB. 
Yet these sources span the typical QSO accretion rates (e.g. 
Steinhardt & Elvis 2010), and indeed our templates are very 
similar to the mean SED in Elvis et al. (1994). It seems that 
these most common QSO SED's are not simply analogous to 
BHB accretion flows. Only the very rare AGN SED's with 
the highest XeiIiI can be well matched to the BHB, as they 
are similar to the high/soft state (see also Done et al. 2011). 



6 DISCUSSION 
6.1 Selection Effects 

Biases and systematics inherent in this sample have already 
been discussed at length in Paper-I and Paper-II. The princi- 
pal imposed selection effect is that sources in our sample are 
bright nearby AGN {z < 0.4). The luminosity of our sources 
is higher than the average among nearby sources, but only 
moderate with respect to samples containing higher redshift 
sources (e.g. GreenOO; LussolO). Our sample contains very 
few sources with L2_iofceV < SxlO**^ erg s~^ or Xsdd < 
0.05. As discussed in Section [3.2.31 and Section [3.6.11 those 
very low luminosity sources may not follow the linear re- 
gression line in the K2-iofceV vs. Xsdd plot in Figure [^ and 
these sources may populate the low L2-iokeV, low K2-iokeV 
region in Figure [TT] It is also possible that these sources may 
not follow other correlations reported in this paper, thereby 
weakening the correlations between L2_iofceV and Mbh and 
H/9 FWHM. Further studies of large samples are required to 
test such possibilities. 

The weak anti-correlation found between Xsdd and 
M_BH (Table [3} also implies some selection effect. Sources 
having both low black hole mass and low mass accretion rate 
are unlikely to be included in our sample as they would be 
too faint. Hence low mass sources in our sample will have 
relatively high Xsdd- For higher mass sources, their lumi- 
nosity will peak when there is a considerable supply of gas 
around them to be accreted. This occurs around redshift 
~2, and as the available gas decreases, the mass accretion 
rate of high mass AGN in the local universe is suppressed 
(so-called downsizing, Fanidakis et al. 2010). Therefore, in 
the local universe high mass sources should have low Xsdd, 
resulting in the weak anti-correlation found between Mbh 
and Xsdd (Done et al. 2011). 

There are also redshift selection effects. Although com- 
parison of parameter correlations with previous work (based 
on larger samples) results in a general consistency, the ques- 
tion remains whether there could be a redshift dependence 
in the correlations we find. However, there are some evi- 
dences that redshift evolution in the spectral properties of 
AGN may not be strong (Fan 2006). 



6.2 Limitations of the Model and Uncertainties 

There is another underlying question, whether the correla- 
tions found might arise artificially as a result of our model 
assumptions. We will consider this point in two ways: 
First, there are no direct constrains on the parameters in 
our SED model. Compared with some previous work (e.g. 
VF07, LussolO, GrupelO), our spectral fitting employs the 
least external constraints on the values of its parameters. 
The only parameter that is directly constrained is Mbh, 
whose value is restricted by the FWHM of the intermediate 
and broad components of the H/3 emission line. However, 
this range often spans more than one order of magnitude, 
and the best-fit Mbh did not exceed the model limits for 
most sources (see Table CI of Paper-I). Therefore this con- 
strain should not be strong enough to cause systematics. 
Second, for previously known correlations such as Xsdd vs. 

K2-10keV, Xsdd VS. r2-10feeV and H/? FWHM vs. r2-10feeV, 

our results are all consistent with past studies based on a 
variety of AGN samples. This suggests that the reported 
correlations should be intrinsic, and that our results are not 
strongly contaminated by model assumptions. As discussed 
in the previous sections, the differences between our results 
and those previously reported are mainly due to two reasons. 
One is our exclusion of highly obscured sources, which re- 
duces the non-intrinsic dispersion within these correlations. 
The other reason is that our parameter values are derived 
from a detailed spectral fitting, rather than from simply ap- 
plying scaling relations which will contain high uncertain- 
ties. 

However, the range of values for each of the 11 key pa- 
rameters discussed previously could be dominated by model 
uncertainties, except for T2-iokeV, whose measurement is 
relatively model independent. For this reason we did not 
adopt the uncertainties returned by Xspec, because they 
must be negligible compared with the model uncertainties. 
Such model uncertainties are very difficult to estimate, and 
the values of the same parameter derived from different SED 
models, can be quite different. As an illustration, in the fol- 
lowing paragraphs we will compare the parameter values be- 
fore and after introducing a colour temperature correction 
into our broadband SED model. 



6.3 The Effect of the Colour Temperature 
Correction 

The colour temperature correction is only important for 
sources having both a low black hole mass and a high mass 
accretion rate (see Done et al. 2011 and references therein). 
So it only affects a small fraction of all the sources in our 
sample, mainly the NLSls. The main consequences of intro- 
ducing a colour temperature correction by using Model-B 
(i.e. optxagnf in Xspec vl2) are that Mbh increases, Lboi 
decreases and so Xsdd decreases. For example, the Xsdd 
of PC 2233-1-134 decreases significantly from 14 to 2.4 af- 
ter using Model-B, making it much less extreme. i-i2-iokeV 
and K5100A also decrease due to the decrease of Ltoi- ctox 
changes slightly but not significantly, because the luminosity 
at 25OOA is mainly constrained by the OM data (Paper-I). 
Figure [T^ compares the distribution of these parameters be- 
tween Model- A (Paper-I) and Model-B fitting. Table [6] lists 
the average values of these parameters for NLSls, BLSls 
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Table 6. The average value of parameters from the best-fit SEDs using Modcl-A (without a colour temperature correction, Paper-I) and 
Model-B (with colour temperature correction). The values and their standard deviations are calculated separately for NLSls, BLSls and 
for the whole sample. 



Model- 


NLSls 
A 


NLSls 
B 


BLSls 
A 


BLSls 
B 


Whole Sample 
A 


Whole Sample 
B 


< K2_10fceV > 


1271 




29lt^ 


30l?? 


41+85 


381^3 


< K5IOOA > 


2Q+37 








16lf 


15l^* 


< ^Edd > 


^•^"^-0.80 




0.27l°:« 


^■^'^-0.17 


nc+1.42 




<Mbh > 


7.11±0.54 


7.37±0.47 


8.04±0.48 


8.10±0.41 


7.83±0.64 


7.93±0.52 


< Olox > 


1.42±0.08 


1.39±0.10 


1.34±0.16 


1.34±0.15 


1.36±0.14 


1.35±0.14 


<Li,o! > 


45.19±0.54 


45.02±0.49 


45.59±0.52 


45.61±0.53 


45.49±0.55 


45.47±0.57 



and the whole sample, for both Model-A and Model-B fit- 
ting. This confirms that the large differences in results from 
the Model-A and Model-B fittings are mostly restricted to 
the NLSls, whose colour temperature corrections are signif- 
icant. 

In order to further investigate the differences between 
using Model-A and Model-B, we redo all of the above cross- 
correlation analysis by adopting Model-A parameter values 
from Paper-I, and then we compare the statistical results 
in the appendices. Appendix [Cl shows all of the correlation 
plots that could in principle be modified by the differences 
between the Model-A and Model-B fittings. In each plot the 
dashed gray lines are for Model-B fitting (this paper), com- 
pared with the solid orange lines for Model-A fitting. It is 
clear that there are no significant changes in any of these 
correlation plots. This is further confirmed by the correla- 
tion matrix for Model-A fitting in Appendix [D] Perform- 
ing a PCA on this matrix, very similar eigenvectors and 
eigenvalues are obtained. Appendix [Elshows the mean SEDs 
based on the parameter values obtained from Model-A fit- 
ting, which does not include the colour temperature correc- 
tion. Therefore, we conclude that use of the refined Model- 
B compared with the original Model-A, does not alter the 
main results reported in this paper, although for individual 
sources such as the NLSls, the refined model should be more 
realistic. 



6.4 The Correction for Radiation Pressure 

Marconi et al. (2008) suggested that the virial mass esti- 
mates should be corrected for the effects of radiation pres- 
sure (RP), especially for those sources with high Eddington 
ratios, such as the NLSls. In our study the M-bh was not de- 
rived directly from the virial mass, but was only constrained 
by the virial masses calculated from the FWHM of the in- 
termediate and broad components of the H/3 line. The final 
estimate of }Abh is derived from the best-fit SED, and so 
there is no need to correct for the RP effect. However, we 
may still choose to derive the virial mass from using the 
FWHM of H/3 profile (narrow component subtracted), and 
then correct these values for the RP effect. Paper-I has listed 
and compared these masses (M_rp) with the best-fit masses 
{Mbf) from the SED fitting using Model-A. The distribu- 
tions of these two estimates of masses are very similar, ex- 



cept that the average M_rp is 0.4 dex higher compared with 
Mbf. 

As discussed in previous section, the SED fittings in- 
cluding the colour temperature correction (Model-B) in- 
creases the average Mbf by 0.27 dex for NLSls, and 0.05 dex 
for BLSls (TableEJ. So for the NLSls using Model-B fitting, 
the average Msf is just 0.05 dex smaller than the average 
Mijp, while for the BLSls the difference is still 0.34 dex. 
Furthermore, no significant difference is found in Table |4] if 
we substitute Mpp for Mpp, and cross-correlate with the 
other 10 SED parameters. The coefficients in Equation 1151 
[16] only differ by less than 1 a when using Mpp instead of 
Mbf- This suggests that the difference between Mpp and 
Mbf is far less than the intrinsic dispersion in any of the 
correlations, and so is not important in our correlation stud- 
ies. 



6.5 The 4000 km s"^ H/? FWHM Break 

The FWHMi^,3 = 2000 km s"^ is the conventional, but 
arbitrary value to distinguish between NLSl and BLSl 
(Goodrich 1989). Recently, the limit of 4000 km s"^ for the 
FWHMh/3 was claimed to be of special interest. For exam- 
ple, when AGNs are divided into two populations based on 
FWHMh/3 = 4000 km s~^ (population A: FWHMh^j < 
4000 km s"^ and B: FWHMff^j > 4000 km s"^), it ap- 
peared that most radio loud sources are contained in popu- 
lation B (Sulentic et al. 2008). Compared with the RL-RQ 
and the NLSl-BLSl divisions, this dividing line in FWHM 
also seems to be more effective in distinguishing the different 
SEDs (Sulentic et al. 2008). Furthermore, ZhoulOa reported 
that in the FWHMj^^ vs. r2_iofcey correlation plot, there is 
a change in the slope at ~4000 km s~^, as consistent with 
our results. They also showed that the the broadest iron Ka 
lines, those with intrinsic width o > 0.5 keV) have are all 
found in AGNs with FWHMj^^ < 4000 km s"\ 

In Section 13.51 we confirmed a slope change at 
FWHMh;3 =i 4000 km s"^ in FWHMh^j vs. V2-wkeV corre- 
lation plot. All the three RL sources in our sample (purple 
symbols in Figure |9l) have FWHMn^j > 4000 km s"^ To 
highlight the 27 population A sources in our sample, we use 
the square symbol to identify these sources in all correlation 
plots, and we show their distribution in the histograms as 
the green region. The 12 NLSls among them are shown as 
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Figure 15. Comparison of parameter distributions between Model-A and Model-B SED fittings. In each panel the shaded cyan histogram 
is based on our modified SED fitting using Model-B (listed in Table lB2t . with the 12 NLSls highlighted by the shaded orange region. 
The solid black line shows the parameter histogram for Model-A fitting (listed in Table 3 of Paper-I) , with the dashed red line indicating 
the 12 NLSls. 



the red region. The two-sided KS-Test was used to deter- 
mine the significance of the difference between the distribu- 
tions of the sub-samples for the NLSl-BLSl division, and 
the population A-B division. Among the 11 key parameters, 
the population A-B provides a somewhat better division for 
the 7 parameters: K2-iofceV, ksiooa, Xsdd, FWHM, aox, 
ax and auv', while NLSl-BLSl division is better for the 
other parameters (see Table [7|). Therefore, the population 
A-B division for our RQ Type 1 AGN sample does not seem 
to be significantly better than the NLSl-BLSl division. 



7 SUMMARY AND FUTURE WORK 

7.1 Summary of Principle Results 

In this third paper in the series, we have studied the SED 
properties of our Type 1 AGN sample. We employ a new 
broadband SED model {optxagnf in Xspec vl2), which in- 
cludes a colour temperature correction, to construct the SED 
for each source in the sample. Various parameters were ob- 
tained from the results of the SED fitting. A detailed sta- 
tistical analysis was performed, which can be divided into 
three major parts: 

• In the first part we studied the diverse correlations 
found among the SED parameters. We divided these param- 
eters into five groups, according to previously reported cor- 
relations. Within each group we conducted a detailed cross- 
correlation analysis, and applied several regression methods. 
Our results are generally in good agreement with previous 



work. However, as a consequence of the unobscured nature 
of the sample and the more reliable parameter values that 
resulted from our refined spectral fitting, we were able to re- 
duce the non-intrinsic dispersion and so obtain the intrinsic 
and better constrained correlations. 

For example, we confirmed the linear correlation in 
Log(L25oo) vs. Log(L2feev) with slope ~1. The correlation 
of Log(«:2_iofcev) vs. Oox can be approximated using a sec- 
ond order polynomial. The correlation between Uox and 
ljOg{XEdd) is weak and dominated by dispersion in intrin- 
sic source properties. Some strong correlations have been 
confirmed e.g. Log(K2-iofcev) vs. Log{XEdd), Log(K5ioo) 

vs. Log(ABdd), Log(K5ioo) vs. Log(K2-10fcev), F2-10keV 
VS. Log(Ai5dd), r2_lOfeeV VS. Log(fi:5loo), Log(Mi3_ff) vs. 

Log(FWHMH;3) and Log(K2-iofeev) vs. Log(MBH). The cor- 
relations in both r2-iofcei' vs. Log(FWHMH/3) and r2-iofceV 
vs. Log(M_Bj/) change slopes as r2_iofce\A decreases to 
-1.8. The break region is FWHMh/j ~ 4000 km s"^ and 
Log(M_B_r/) ~ 8.0. r2_iofeeV is almost independent of the 
FWHMh/3 after the break region, with a mean value of 1.8. 
By presenting the correlations between ax, auv, Oox and 
L2~iokeV, we also confirmed that AGNs with bluer opti- 
cal/UV spectra tend to have steeper X-ray spectra, and 
their hard X-ray emission is also weaker. Other strong (anti- 
)correlations were found for ax and auv, such as ax (ceuv) 
vs. K2-iokeV, ax (auv) vs. K5100A, ax (auv) vs. FWHMh/3 
and ax {auv) vs. Mbh- 

• In the second part of our work, we performed a sys- 
tematic cross-correlation study by producing the correla- 
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Table 7. The significance level of the difl'ercnce in the NLSl-BLSl pair (Nl-Bl) and Population A-B pair (Pop. A-B) for values of the 
11 key parameters. The two-sided Kolmogorov-Smirnov test is applied. A smaller value suggests a greater difference within each sample 
division pair. 
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tion matrix of the 11 key parameters: r2-iokeV, H2-iokeV, 
kbiooa, Afidd, FWHMfr/3, Mbh, ctox, Lto!, L2-iofceV, ax and 
auv- The PC A was performed on the correlation matrix 
to discover the principal eigenvectors that drive the most 
correlations. Wc found that the first two eigenvectors (PCI 
and PC2) contain ~80% of all correlations in the matrix. 
PCI strongly correlates with Mbh, while PC2 is dominated 
by Ltoi- In addition both PCI and PC2 well correlate with 
^Edd- Thus the two principle eigenvectors are driven by 
Mbh, \Edd and Li,oi (or equivalently M). Our eigenvectors 
also have similar properties to the two principal eigenvectors 
derived by Boroson (2002) based on their optical emission 
line study. 

• In the third part we produced various mean SEDs 
classified by each of the key parameters. The SED shapes 
are found to exhibit similar changes with most parame- 
ters except Lfcoi • This explains the strong correlations found 
among these key parameters. A more detailed mean SED 
comparison suggests that the dispersion within each of the 
three binned SEDs is minimized (and the difference between 
them is maximized) for composite parameters which depend 
on both \Edd and Mbh, such as K2-wkeV, ksioo and 
FWHM. This is because the SED change is not determined 
solely by any one of these key parameters. It should ulti- 
mately depend on both Asdd and Mbh- 

• To test the robustness of the main results from our 
correlation study, we looked at the black hole masses cor- 
rected for the effect of radiation pressure. We found no sig- 
nificant differences from using our best-fit black hole masses. 
We also compared the correlation results between Model-A 
(without a colour temperature correction) and Model-B (in- 
cluding a colour temperature correction) fitting, and found 
that they were all very similar. 

• The population A-B division for AGNs was compared 
with the NLSl-BLSl division, but it did not prove to be a 
better AGN classification method. 



7.2 Future Work 

Our sample is limited to relatively high A^dd, with few ob- 
jects below Afidd = 0.05. These (predominantly LINER) ob- 
jects are the ones expected to bo the counterparts of the 
low/hard state in BHB. Another important extension would 
be to increase the sample size and include rare higher mass 
objects with high Eddington ratios. This would allow the 
SEDs to be co-added for different Asdd at a given (fixed) 
black hole mass, thus providing a direct comparison with 
the BHB states. 

The major result of this study is that the SEDs of AGN 
exhibit a very wide range, most plausibly as a function of 
mass accretion rate for a given mass black hole. This clearly 
shows that so-called unified schemes, where AGN have in- 



trinsically identical spectra which are modified by orienta- 
tion dependent obscuration, are an over simplification of the 
actual situation. In fact, unobscured AGN can have quite 
different SED shapes depending on Xsdd, and Mbh- 

Although not widely appreciated, this is broadly ex- 
pected by analogy of AGN with BHB. The stellar mass black 
holes clearly show a dramatic change in spectral shape with 
XEdd, but unlike AGN, these changes can be tracked in a 
single object because of the much shorter timescale for vari- 
ability. However, while the highest Xsdd spectra appear sim- 
ilar to the disc dominated 'high/soft state' seen in BHB, the 
more typical AGN (with an SED similar to the standard 
quasar SED template in Elvis et al. (1994)), do not appear 
to have SED properties which match with any spectral state 
known in BBB. 

This might indicate a physical break in the properties 
of the accretion flow between stellar mass and supermassive 
black holes. The most obvious change in physical conditions 
between these two mass regimes is that AGN discs arc strong 
in the UV, and so are capable of powering substantial mass 
loss via a UV line driven wind. A consequence of mass loss 
in the wind is that the accretion rate is no longer constant 
as a function of radius, causing an intrinsic change in the 
structure of the accretion flow (e.g. Proga, Stone & Kali- 
man 2000). Emission/absorption/scattering processes in the 
wind can also change the observed properties of the spec- 
trum (Sim et al. 2010). Further work on theoretical disc 
models including these effects will show whether standard 
AGN accretion flows are indeed sculpted by a wind. 
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Table Al. A summary of references about the correlations among the most important AGN SED parameters. The upper right triangle 
region shows some most recent works about each correlations pair, while the lower left triangle region shows which cross-correlations 
are studied in this paper, 'a/' means this cross-correlation pair is studied in this paper, 'x' means not studied. The 'FWHM' is the 
FWHM of narrow component subtracted H/3 profile. 
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(RM)*: the reverberation mapping studies, e.g. Kaspi et al. (2000), Peterson et al. (2004); GreenOQ: Green et al. (2009); G99: Grupe 
et al. (1999); GIO: Grupe et al. (2010); H07: Hopkins, Richards & Hcrnquist (2007); FIO: Fanidakis et al. (2010); Paper-I: Jin et al. 
2011; LIO: Lusso et al. (2010); M04: Marconi et al. (2004); S06: Shemmer et al. (2006); SOS: Shemmer et al. (2008); V07: Vasudevan 
& Fabian (2007); V09: Vasudevan & Fabian (2009); Woo02: Woo & Urry (2002); ZhoulOa: Zhou & Zhang (2010); ZhoulOb: Zhou & 
Zhao (2010). 
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APPENDIX B: PARAMETER VALUES FROM MODEL-B FITTING {OPTXAGNF: WITH COLOUR 

CORRECTION) 
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Table Bl. Broadband SED fitting parameters using Modcl-B, and the fitting outputs {Lbol^ fd^ fc, fp)- ID: object number, 

the same as Table 1 in Papcr-I; T^u.gal Njj.int: the fixed galactic and free intrinsic neutral hydrogen column densities 

in 10'^" cm^'^; Tpow'- the powcrlaw component's slope in the SED fitting, (*) denotes the objects whose powerlaw slopes 

hit the uplimit of 2.2 and were fixed there; Fpl: the fraction of powcrlaw component in the total reprocessed disc emission; 

Ticor- corona (truncation) radius in unit of Gravitational radii (ig) within which all disc emission is reprocessed into the 

Comptonisation and powerlaw components; Te- temperature of the Compton up-scattering electron population; Tau: optical 

depth of the Comptonisation component; log{Msjj): the best-fit black hole mass; log(M): total mass accretion rate; Lj^;: 

bolometric luminosity integrated from 0.001 keV to 100 keV; fd, fc, fp- luminosity fractions of disc emission, soft Comptonisation 

and hard X-ray Compotonisation components in the bolometric luminosity; x^' the reduced of the broadband SED fitting. 

,int ^pow Fpl Rcor Te Tau logiMBH) lo9{M) fa fc fp 

xlO^° xlO^° Tg keV Mq g s-1 JO"*^ reduced 

100. 0.246 17.4 8.61 26.08 62.2 0.19 0.25 0.56 1.00 

100. 0.212 16.3 7.84 25.25 9.11 0.19 0.49 0.32 0.97 

10.2 0.214 12.2 7.61 25.92 42.8 0.87 0.10 0.03 1.14 

100. 0.317 15.2 8.78 25.45 14.4 0.19 0.41 0.40 1.16 

69.0 0.202 14.4 7.94 26.37 119 0.32 0.44 0.24 1.04 

20.6 0.348 12.2 8.50 26.41 130 0.59 0.22 0.19 1.09 
10.5 0.149 33.4 7.58 26.01 52.9 0.26 0.53 0.21 1.16 

23.7 0.598 6.91 7.54 25.44 14.3 0.58 0.35 0.06 1.42 
16.5 0.267 15.0 7.40 25.70 25.9 0.84 0.05 0.11 20.7 
100. 0.295 13.7 8.34 26.24 89.9 0.19 0.55 0.26 1.40 
22.5 0.775 8.40 7.92 26.02 54.2 0.65 0.18 0.17 0.96 
19.5 0.385 11.0 7.82 25.32 10.7 0.65 0.18 0.17 1.13 
11.4 0.142 17.9 9.20 26.50 161 0.90 0.01 0.09 1.33 
100. 0.294 16.0 8.24 25.85 36.0 0.19 0.23 0.58 1.07 
100. 0.764 4.33 6.30 25.28 9.72 0.80 0.18 0.02 3.67 
100. 0.242 15.1 8.30 25.93 43.9 0.19 0.56 0.25 1.40 
88.9 0.474 8.27 7.79 26.12 68.4 0.35 0.56 0.09 1.06 
11.9 0.260 13.6 8.12 26.33 110 0.51 0.37 0.12 1.48 

31.4 0.144 20.6 7.71 24.92 4.28 0.37 0.52 0.12 1.03 

12.1 0.186 22.7 6.80 24.87 3.81 0.94 0.04 0.02 1.55 
23.1 0.205 19.5 7.98 26.09 63.0 0.61 0.10 0.29 3.45 

47.7 0.116 22.1 7.84 25.50 16.2 0.34 0.52 0.14 1.09 

95.8 0.628 9.70 8.00 25.06 5.92 0.22 0.43 0.35 0.99 
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Table B2. Key Parameters from Modcl-B fitting. ID: object number, the same as Table 1 in Paper-I; T2—u)keV- the slope of the single 



powerlaw fitted to 2-fO keV speetrum. 1^2—iOkeV- 2-10 keV luminosity (in 10 erg s 



coefficient; AL„ 



the monochromatic luminosity at 2500A (in IC* 



erg s 



i^2-u)keV- the 2-lOkeV bolometric correction 
/: the monochromatic luminosity at 2keV (in lO**^ 



erg s~^); aox' the optical X-ray spectral index; ALsioq; the monochromatic luminosity at 5100A (in 10** erg s~^); k^ioo- the 5100A 
bolometric correction coefficient; FWHMjj^: the narrow component subtracted H/3 FWHM; l^bol/^Edd- the Eddington Ratio; ax- the 
soft X-ray slope between 0.2-2 keV (corrected for Galactic and intrinsic absorption), assuming Fi, oc i/~"; auv- the optical/UV slope 
between 1700-6500A (corrected for Galactic and intrinsic reddening and absorption). 
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APPENDIX C: PARAMETER CROSS-CORRELATIONS USING VALUES FROM MODEL-A FITTING 

{OPTXAGN: WITHOUT COLOR CORRECTION) 

Similar cross-correlation plots as reported in previous sections but use parameter values from Model-A fitting in Paper-I. In 

each panel, the various point symbols show Model-A data. The solid orange line is the regression result for these Model-A 
data using the same regression methods as in previous sections. The dashed gray line is our result from previous sections 
based on Model-B data. All dotted lines are from the other literatures as indicated in previous sections. 
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Table Dl. The cross-correlation matrix of the 9 key parameters based on Model-A data from Jin et al. (2011). ID 1~9 are given to 
each parameter, ps is the Spearman's rank coefficient, is probability of random distribution in logarithm, a and /3 are the bisector 
regression line coefficients assuming Y=/3X+a. The coefficients in the upper right triangle region assumes the vertical parameters to 
be X. The coefRcients in the lower left triangle region assumes the horizontal parameters to be X. 
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APPENDIX E: MEAN AGN SEDS USING MODEL-A FITTING (OPTXAGN: WITHOUT COLOR 

CORRECTION) 
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Figure El. The AGN mean SEDs based on different values of the 9 key parameters from Model- A fitting in Paper-I (i.e. without effect 
of color temperature correction). For each parameter, the 51 sources are sorted according to the parameter value, and are then divided 
into three subsets evenly so that each subset contains 17 sources. PG 1004+130 is excluded from its subset. Finally, a mean SED is 
constructed for each of the three subsets after renormalizing each individual SED to the mean 2500A luminosity of that subset. The 
three panels (A, B, C) in each row show the mean SEDs of the subsets classified by the parameter shown in the panel title. In each 
panel, the solid curve is the mean SED, while the color shaded region is the ztltr deviation. The 2500 A and 2 keV positions are marked 
by the vertical solid orange lines, whose related height shows the value of aox, The SED peaking position is also marked by the vertical 
solid purple line. The average values of some other parameters in that subset are also shown in the panel. All the mean SEDs have been 
rescaled by the same arbitrary constant in the Y-axis which is 1.3xl0~*®. Note that the energy ranges 1.4 eV < E < 6 eV and 0.3 keV 
< E < 10 keV are covered by SDSS, OM and EPIC data, while the SED in the rest energy bands is determined by the extrapolating of 
the best-fit model. 
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Figure El. continued 
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